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Abstract 

The interest in X/7-ray Astronomy has grown enor- 
mously in the last decades thanks to the ability to send 
X-ray space missions above the Earth's atmosphere. 
There are more than half a million X-ray sources de- 
tected and over a hundred missions (past and currently 
operational) devoted to the study of cosmic X/7 rays. 
With the improved sensibilities of the currently active 
missions new detections occur almost on a daily ba- 
sis. Among these, neutron-star X-ray binaries form an 
important group because they are among the brightest 
extra-solar objects in the sky and are characterized by 
dramatic variability in brightness on timescales rang- 
ing from milliseconds to months and years. Their main 
source of power is the gravitational energy released by 
matter accreted from a companion star and falling onto 
the neutron star in a relatively close binary system. 

Neutron-star X-ray binaries divide into high-mass 
and low-mass systems according to whether the mass 
of the donor star is above ^8 or below ^2 Mq, respec- 
tively. Massive X-ray binaries divide further into super- 
giant X-ray binaries and Be/X-ray binaries depending 
on the evolutionary status of the optical companion. 
Virtually all Be/X-ray binaries show X-ray pulsations. 
Therefore, these systems can be used as unique natu- 
ral laboratories to investigate the properties of matter 
under extreme conditions of gravity and magnetic field. 

The purpose of this work is to review the obser- 
vational properties of Be/X-ray binaries. The open 
questions in Be/X-ray binaries include those related to 
the Be star companion, that is, the so-called "Be phe- 
nomenon" , such as, timescales associated to the for- 
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mation and dissipation of the equatorial disc, mass- 
ejection mechanisms, V/R variability, and rotation 
rates; those related to the neutron star, such as, mass 
determination, accretion physics, and spin period evo- 
lution; but also, those that result from the interaction 
of the two constituents, such as, disc truncation and 
mass transfer. Until recently, it was thought that the 
Be stars' disc was not significantly affected by the neu- 
tron star. In this review, I present the observational 
evidence accumulated in recent years on the interaction 
between the circumstellar disc and the compact com- 
panion. The most obvious effect is the tidal truncation 
of the disc. As a result, the equatorial discs in Be/X- 
ray binaries are smaller and denser than those around 
isolated Be stars. 

Keywords X-rays: binaries - stars: neutron - stars: 
binaries close -stars: emission line. Be 

1 Definition and classification of X-ray binaries 

In very general terms, one can simply define X-ray bina- 
ries as systems that consist of a compact object orbiting 
an optical companion. They are " close" binary systems 
because there exists a transfer of mass from the optical 
component to the compact object. By "optical compan- 
ion" it is understood that nuclear burning is still taking 
place in its interior. Figure [T] shows a tree-diagram de- 
picting all the different subsystems that comprise the 
generic group of X-ray binaries. 

In referring to the two components in an X-ray bi- 
nary one should be careful and learn which is the sub- 
ject of investigation as the same name can be used to 
mean different components. In massive X-ray bina- 
ries, the most massive star is normally termed primary 
whereas the less massive one is called secondary. In low- 
mass systems, the term primary refers to the neutron 
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Fig. 1 Classification of X-ray binaries. 

star while the word secondary is reserved for the late- 
type companion. Other names also used are "optical" 
or "donor" for the larger star and "compact", "gainer" 
or "accreting" for the denser companion. 

They admit several classification schemes depend- 
ing upon whether the emphasis is put on the type of 
the compact companion or the physical properties of 
the optical star. X-ray binaries divide up into black- 
hole systems, neutron star X-ray binaries or cataclysmic 
variables (if the compact object is a white dwarf) . Nev- 
ertheless, the term "X-ray binaries" is normally re- 
served to designate binaries with neutron stars. 

1.1 High- mass X-ray binaries 

Neutron-star X-ray binaries are divided up into high- 
mass (HMXBs) and low-mass (LMXBs) X-ray binaries 
depending on the spectral type of the mass donor, as 
this feature determines the mode of transferring mass 
to the compact object and the environment surround- 
ing the X-ray source. HMXBs contain early-type (O or 
B) companions, while the spectral type of the optical 
star in LMXBs is later than A. HMXBs are strong emit- 
ters of X-ray radiation. Sometimes they appear as the 
brightest objects of the X-ray sky. The high-energy ra- 
diation is produced as the result of accretion of matter 
from the optical companion onto the neutron star. The 
term accretion refers to the gradual accumulation or de- 
position of matter onto the surface of an object under 
the influence of gravity. If the accreting object is a neu- 
tron star (or black hole), then matter falls down onto 
an enormous well of gravitational potential and is accel- 
erated to extremely high velocities. When the matter 
reaches the surface of the neutron star, it is rapidly de- 
celerated and the free-fall kinetic energy radiated away 
as heat which is available to power the X-ray source. 

The luminosity class serves to subdivide HMXBs into 
Be/X-ray binaries (BeXB), when the optical star is a 
dwarf, subgiant or giant OBe star (luminosity class HI, 



IV or V) and supergiant X-ray binaries (SGXBs), if 
they contain a luminosity class I-II star. 

In SGXBs, the optical star emits a substantial stel- 
lar wind, removing between 10"^ - 10"^ Mq yr'^ with 
a terminal velocity up to 2000 km s^^. A neutron star 
in a relatively close orbit will capture some fraction of 
this wind, sufficient to power a bright X-ray source. 
If mass transfer occurs via Roche lobe overflow, then 
the X-ray emission is highly enhanced and an accretion 
disc is formed around the neutron star. At present, 
there is known only one disc-fed SGXB in the Galaxy 
(Gen X-3) and three in total (SMG X-1 and LMG X-4), 
while there are about a few tens of wind-fed SGXBs. 
Because of their brightness and persistent X-ray emis- 
sion, SGXBs were the first to be discovered. They were 
initially thought to represent the dominating popula- 
tion of HMXBs, whereas BeXB were considered atyp- 
ical cases. Hence, the name classical or standard was 
given to SGXBs. 

In BeXB, the optical companion is a Be star. Be 
stars are non-supergiant fast-rotating B-type and lu- 
minosity class III-V stars which at some point of their 
lives have shown spectral lines in e mission, hence the 
qualifier "e" in thei r spectral types ( Porter fc Riviniu^ 
20031: iBalonal [200o[ Isiettebaki flissl ) . The best stud- 



ied lines are those of hydrogen (Balmer and Paschen 
series) but Be stars ma.v also show He, Fe in emission 
(see e.g Hanuschi^ll996l ). They also show an amount of 
IR radiation than is larger than that expected from an 
absorption-line B star of the same spectral type. This 
extra long- wavelength emission is known as infrared ex- 
cess. 

The origin of the emission lines and infrared excess 
in BeXB is attributed to an equatorial disc, fed from 
material expelled from the rapidly r otating Be star in a 
mann er that it is not yet understood ( Porter fc Riviniui 
20031) . During periastron, the neutron star passes close 



to this disc, sometimes may even go through it caus- 
ing major disruption. A large flow of matter is then 
accreted onto the compact object. The conversion of 
the kinetic energy of the in- falling matter into radiation 
powers the X rays. BeXB have large orbital periods and 
by definition they are non-supergiant systems. Hence, 
the Be star is well within the Roche lobe. However, 
transient Roche lobe overflow may occur during perias- 
tron passage in systems with large eccentric orbits or 
during giant X-ray outbursts when a large fraction of 
the Be star's disc is believed to be accreted. 

Most BeXB are transient systems and present mod- 
erately eccentric or bits (e > 0.3), althou gh persistent 
sources also exist ( Reig fc Roche 19991) . Persistent 
BeXB differ from transient BeXB in that they display 
much less X-ray variability (no large outbursts are de- 
tected), lower luminosity (L^ < 10'^^ erg s""'^), contain 
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Table 1 List of galactic BcXB. Only systems with known optical counterparts and well-established optical/X-ray behaviour are included 



X-ray name 


Optical Ctp. 


RA (J2000) 


Dec (J2000) 


Spec, type 


V 


J 


E(B-V) 




P„b (d) 


e 


d (kpc) 


4U 0115-f-634 


V635 Gas 


01 


18 


31 


90 


+63 44 24.0 


BO. 2 Vo 


14.8-15.5 


10.8-12.3 


1.55 


3.6 


24.3 


0.34 


8 


IGR J01363+6610 




01 


36 


18 


00 


+66 10 36.0 


Bl IV- Ve 


13.3 


10 


1.61 








2 


RX J0146. 9-1-6121 


LS I +61 235 


01 


47 


00 


17 


+61 21 23.7 


BlVo 


11.2 


10 


0.93 


1412 


330 




2.3 


IGR J01583-I-6713 




01 


58 


18 


20 


+67 13 26.0 


B2IVo 


14.4 


11.5 


1.46 


19692 






4 


RX J0240. 4+6112 


LS I +61 303 


02 


40 


31 


67 


+61 13 45.6 


BO. 5 Vo 


10.7 


8.8 


0.75 




26.45 


0.54 


3.1 


V 0331-1-530 


BQ Cam 


03 


34 


59 


89 


+53 10 23.6 


08-9 Vo 


15.1-15.8 


11.4-12.2 


1.9 


4.4 


34.3 


0.3 


7 


4U 0352-1-309 (X For) 


HD 24534 


03 


55 


23 


08 


+31 02 45.0 


O9.5IIIc-B0Vo 


6.1-6.8 


5.7-6.5 


0.4 


837 


250 


0.11 


1.3 


RX J0440. 9-1-4431 


LS V +44 17 


04 


40 


59 


32 


+44 31 49.3 


B0.2Vo 


10.8 


9.2 


0.65 


203 






3.3 


lA 0535-1-262 


V725 Tau 


05 


38 


54 


57 


+26 18 56.8 


BO IIIo 


8.9-9.6 


7.7-8.5 


0.75 


105 


111 


0.47 


2.4 


IGR J06074+2205 




06 


07 


24 


00 


+22 05 00.0 


B0.5V 


12.3 


10.5 


1.88 








5 


XTE J0658-073 




06 


58 


42 


00 


-07 11 00.0 


09.7 Ve 


12.1 


9.7 


1.19 


160.4 






3.9 


4U 0726-260 


V441 Pup 


07 


28 


53 


60 


-26 06 29.0 


08-9Ve 


11.6 


10.4 


0.73 


103.2 


34.5 




6 


RX J0812. 4-3114 


LS 992 


08 


12 


28 


84 


-31 14 52.2 


BO. 2 Ill-IVe 


12.4 


11.2-12.0 


0.65 


31.89 


80 




9 


GS 0834-430 




08 


35 


55 


00 


-43 11 06.0 


BO-2 III-Vo 


20.4 


13.3 


4.0 


12.3 


105.8 


0.12 


5 


GRO J 1008-57 




10 


09 


46 


00 


-58 17 30.0 


09e-Blo 


15.3 


10.9 


1.9 


93.5 


247.5 


0.66 


2 


RX J1037. 5-5647 


LS 1698 


10 


37 


35 


50 


-56 48 11.0 


BO III-Vo 


11.3 




0.75 


862 






5 


lA 1118-615 


Hen 3-640 


11 


20 


57 


21 


-61 55 00.3 


09.5 Ill-Ve 


12.1 


9.6 


0.92 


405 






5 


4U 1145-619 


V801 Con 


11 


48 


00 


02 


-62 12 24.9 


Bl Vno 


9.3 


8.7 


0.29 


292.4 


187.5 


>0.5 


3.1 


4U 1258-613 


GX 304-1 


13 


01 


17 


20 


-61 36 07.0 


B2 Vno 


13.5 


9.8 


2.0 


272 


132.5 


>0.5 


2.4 


2S 1417-624 




14 


21 


12 


80 


-62 41 54.0 


Bl Vo 


17.2 


13.3 


2.0 


17.6 


42.12 


0.45 


10 


GS 1843+00 




18 


45 


36 


90 


+00 51 48.2 


BO-2 IV- Ve 


20.9 


13.7 


2.5 


29.5 






>10 


XTE J 1946+274 




19 


45 


39 


30 


+27 21 55.0 


BO-1 IV-Ve 


16.8 


12.5 


2.0 


15.8 


169.2 


0.33 


8-10 


KS 1947+300 




19 


49 


30 


50 


+30 12 24.0 


BO Ve 


14.2 


11.7 


1.09 


18.76 


40.4 


0.03 


9.5 


EXO 2030+375 


V2246 Cyg 


20 


32 


15 


20 


+37 38 15.0 


BOe 


19.7 


12.1 


3.8 


41.8 


46.03 


0.41 


5 


GRO J2058+42 




20 


59 


00 


00 


+41 43 00.0 


O9.5-B0 IV- Vo 


14.9 


11.7 


1.38 


192 


110 




9 


SAX J2103. 5+4545 




21 


03 


35 


71 


+45 45 05.5 


BO Vo 


14.2 


11.8 


1.35 


358.6 


12.67 


0.4 


6.8 


4U 2135+57 


Cep X-4 


21 


39 


30 


72 


+56 59 10.0 


B1-B2 Vo 


14.2 


11.8 


1.3 


66.3 






3.8 


SAX J2239. 3+6116 




22 


39 


20 


90 


+61 16 26.8 


BO-2 III-Vo 


15.1 


11.5 


1.4 


1247 


262.6 




4.4 
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Table 2 Statistics on HMXBs in the Milky Way 

Number of neutron-star X-ray binaries^ 327 

Number of suspected HMXE^ 131 

Number of suspected BeXB^ 63 

Number of confirmed BeXB'^ 28 



^Sources in tIie | Liu et al. and lLiu et al. 1 1 I2OO6I) catalogs plus update from the 4th IBIS/ISGRI soft gamma-ray survey catalog 

l|Bird et al.ll201of i 

Sources in the updated on-line version of the iRaguzova fc Popovl l l2005h catalog (Iht tp: / /xray.sai.msu.ru/~raguzova/BeXcat/[ | 
^ Systems with known optical counterpart and verified X-ray behaviour (from Table 



of Corbet's diagram. Disc-fed SGXBs (squares) show 
short orbital periods and short spin periods and dis- 
play an anticorrelation in the Porb — -Pspin diagram. The 
small orbital separation and evolved companions make 
Roche lobe overflow the most likely mass transfer mech- 
anism. Wind-fed SGXBs (triangles) show long spin pe- 
riods and short orbital periods, occupying a more or 
less flat region in the Porb — Pspin diagram. Two sys- 
tems (open triangles) prevent the region from being a 
horizontal line: OA01657-41, which might be making 
a transition to the disc-fed SGXB and 2S 0114-1-65, for 
which the association of the lO'' s pulsations with the 
spin period of the neutron star remains controversial 
( Koenigsberger et al.ll2006l ). The spin and orbital peri- 
ods of BeXB (filled circles) exhibit a clear correlation. 
The open circle in Fig. [2] represents SAX J2103. 5-1-4545, 
a peculiar system whose X-ray properties are similar to 
those of wind-fed systems b ut whose optical /I R emis- 



sion resemble that of BcXB (jReig et al. 



l2010ah . 



of a neutron star and a Be star. However, there seems 
to be no apparent mechanism that would prevent the 
formation of Be stars with black holes (BH) or white 
dwarfs (WD). Surprisingly, not a single BeXB is known 
to host a black-hole companion in our Galaxy, whereas 
the interpretation of 7 Gas as a Be-(-WD system still 
remains very controversial. 

Two ideas have been put forw ard to expla i n this 
apparent lack of Be/BH binaries. IZhang et al.l (j2004[ ) 
extended the application of the v is cous decretion disc 
model of lOkazaki fc Negueruelal (j200l[ ) to compact 
companions of arbitrary mass and showed that the most 
effective Be disc truncation would occur in relatively 
n arrow systems. Using the po pulation synthesis results 
of iPodsiadlowski et al.l (|2003( ). which state that binary 
black holes are most likely to be born in systems with 
narrow orbits (Porb ^30 days), the reason for the lack of 
these systems can be attributed to the difficulty to de- 
tecting them. Be/BH binaries are expected to be X-ray 
transients with very long quies cent states. In contrast, 
Belczvnski fc ZiolkowsM |2009t ) showed that Be/BH bi- 



1973t IStefla et al.lll986t IWaters fc van Kerkwiiklll989l) . 



The observed correlation in Be/X systems is ex- 
plained in terms of the equilibrium period, defined as 
the period at which the outer edge of t he magneto- narie s do not necessarily have narrow orbits. These 
sphere rotates with the Keplerian velocity (Davidson fc Ostriktet^ rs argued that the predicted ratio of Be/NS bina- 
ries to Be/BH binaries (FNS-to-BH ~ 30 — 50) in our 
Galaxy, based on current population synthesis models 
and evolutionary scenarios, is consistent with the obser- 
vations: there are 60 known Be/NS binaries, hence one 
would expect 0-2 Be/BH bina ries, consistent with the 
observed galactic sample . BothlZhang et al. ( 2004 ) and 
Belczvnski fc Ziolkowskil |2009l ) agree in that Be/BH 



If the the neutron star (and hence the magnetosphere) 
rotates faster that the equilibrium period, then matter 
is spun away by the propeller mechanism; only when 
the spin period is larger than the equilibrium period can 
matter be accreted on to the neutron star surface. This 
results in angular momentum transfer to the neutron 
star, increasing its rotation velocity (decreasing the spin 
period). The equilibrium period depends mainly on the 
mass flux (or accretion rate) because it determines the 
size of the magnetosphere which is assumed to corotate 
with the neutron star. In turn, the accretion rate de- 
pends on the separation of the two components of the 
binary systems, hence on the orbital period. 

The compact object in all confirmed BcXB (Ta- 
ble [OJ is a neutron star. In fact, many times the neu- 
tron star is taken as a defining property of BeXB. It is 
common to read in the literature that a BeXB consists 



binaries should exist in the Galaxy. 

1.2 Other types of high-mass X-ray binaries 

The traditional picture of two classes of HMXBs, 
namely SGXB — subdivided into low-luminosity or 
wind-fed systems and high-luminosity or disc-fed sys- 
tems — and BeXB — transient and persistent — is 
giving way to a more complex situation where newly 
discovered systems may not fit in these categories. 
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1.2.1 Low eccentricity BeXB 

There is a group of so far five BeXB (X Per, GS 0834- 
430, KS 1947+300, XTE J1543-568, and 2S 1553-54^1), 
characterised by Poib^ 30 d and very low eccentricity 
(e < 0.2). Their low eccentricity requires that the neu- 
tron star received a much lower kick velocity at birth 
tha n previously assum ed by current evolutionary mod- 
els |Pfahl et al.ll2002[ ). 

Most popular models for neutron star kicks involve 
a momentum impulse delivered around the time of the 
core collapse that produced the neutron star. These 
models assume that neutron stars are born with speeds 
in excess of 100-200 km s^^. Such velocities imply a 
probability close to zero that the post-supernova ec- 
centricity is less than 0.2. Hence, the discovery of 
these low-eccentricity BeXB with wide orbits was un- 
expected. Tidal circularizacion is ruled out since this 
mechanism requires that the star almost fills its Roche 
lobe and that there be an effective mechanism for damp- 
ing the tide. Tidal torques should have little effect on 
the orbit of a HMXB with Porb = 10 days, as long as 
the secondary is not too evolved and the eccentricity 
is not so large that the tidal interaction is enhanced 
dramatically at periastron (as it is the case of typical 
B eXB). 



Pfahl et al. 



(120021 ) developed a phenomenological 
model that simultaneously accounts for the long-period 
(>30 d), low-eccentricity (<0.2) HMXBs and which 
does not violate any previous notions regarding the 
numbers and kinematics of other neutron star popu- 
lations. They propose that a neutron star receives a 
relatively small kick (^50 km s~^) if the progenitor, i.e., 
the core of the initially more massive star in the binary, 
rotates rapidly. This condition may be met when the 
progenitor star experienced case Be or Ce mass trans- 
fer in a binary systerrH. If the hydrogen-exhausted core 
of an initially rapidly rotating massive star is exposed 
following case Be or Ce mass transfer in a binary, then 
the core is also likely to be a rapid rotator. 

1.2.2 Obscured sources & Supergiant fast X-ray 
transients 

Since the launch of INTEGRAL in October 2002, the 
situation is also changing among the SGXB group. 



iThe optical counterparts of XTE J1543-568 and 2S 1553-542 
are not known, hence they do not appear in Tabic [TTT1 

^In case B, mass transfer occurs during the shell hydrogen- 
burning phase, but prior to central helium ignition, while case 
C evolution begins after helium has been depleted in the core. 
Cases B and C are naturally divided into an early case (B^ or 
Ce), where the envelope of the primary is mostly radiative, and 
a late case (B; or Cj), where the primary has a deep conveetive 
envelope. 




Fig. 3 Spectral energy d istribution of LSI -1-61 303. From 
IChernvakova et al.l (|2006l ') 



INTEGRAL has unveiled a population of highly ob- 
scured HMXBs with supergiant companions and a new 
type of source displaying outbursts which are signif- 
icantly shorter than typical for BeXB and which are 
characterized by bright flares with a duration of a 
few hours and peak luminosities of lO'^^ — 10^^ erg 
s^^. These new systems have been termed as Super- 



Negucruela et al. 


2006 


Neeucruela et al. 


20081) 



giant fast X-ray Transients (SFXTs. ISmith et al. 

200d: IWalter fc Zurita-Heraj 



2006 



2007 



Both obscured HMXBs and SFXTs display X-ray 
and IR spectra typical of SGXBs. In some cases, 
the X-ray sources are pulsed and orbital parameters 
typical of persistent SGXBs have been found (e.g., 
Bodaghee et~all l2006t [Zurita Heras et al.l l2006f ). The 
heavily-absorbed sources had not been detected by pre- 
vious missions due to high absorption, which renders 
their spectra very hard. The current understanding is 
that the entire binary system is surrounded by a dense 
and absorbing circumstcllar materi al envelope o r co- 
coon, made of cold gas and/or dust (jChatvl 120081 ) . 

SFXTs differ from SGXBs because they are only 
detected sporadically , during very brief outbursts 
( Romano et al. 201ol ). A promising model to ex- 



plain SFXTs invokes highly structured (clumpy) stellar 
winds. The outburst occurs as a result of the accretion 
of one of the clumps of dense matter from the wind. An 
alternative model ("Bc-type" model) assumes a very el- 
liptical orbit for the binary. In this model outbursts are 
triggered when the compact object travels through its 
periastron. Other possibilities imply that SFXTs con- 
tain strongly magnetised neutron stars. The outbursts 
result from t he overcoming o f centrifugal and magnetic 
barriers (see Grebenev 2010l and references therein). 



1.2.3 "f-ray binaries 

7-ray binaries are HMXBs that emit most of their ra- 
diative output in the MeV-TeV range. Currently, only 
four HMXBs are well-established members of this group 
of high-energy sources: LS I +61 303, LS 5039, PSR 
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B1259-63, and Cygnus X-3, while other two are firm 
candidates: Cygnus X-1, and HESS J0632+05711. The 
optical counterpart is either a luminous 0-type star 
(Cygnus X-1, LS 5039), a Be star (LS I +61 303, PSR 
B1259-63, HESS J0632+057) or a likely Wolf-Rayet 
star in Cygnus X-3. In addition to the variable 7-ray 
emission these systems share in common a resolved ra- 
dio counterpart with a jet or jet-like structure, multi- 
wavelength orbital modulation and spectral energy dis- 
tribution (Fig. [3l). The wide range of orbital parameters 
( Paredes 20081 ) and the non- unique nature of the com- 
pact companion (unknown in most systems but with a 
confirmed neutron star in PSR B1259-63 and a fairly se- 
cure black hole in Cyg X-1) represent a challenge for the 
theoretical modelling of these systems. Two alternative 
scenarios may explain the variable 7-ray emission: the 
microquasar or accretion-powered scenario and the pul- 
sar wind scenario. 

All confirmed 7-ray binaries show a jet or jet-like 
radio structure, which would indicate the presence of 
relativistic particles. In black-hole binaries the radio 
jet can account for their broad-ban d spectrum, from 
radio to X-rays (jMarkoff et al.ll2003l) as well as for the 



origin of most of the timing variability (jKvlafis et al, 
2008[ ). Therefore, it is reasonable to think that the 



jet may also be the origin of the very high-energy 7- 
rays. In the context of relativistic jets, the most efficient 
gamma-ray mechanism would be inverse Compton scat- 
tering, by which relativistic particles collide with stellar 
and/or synchrotron photons and boost their energies to 
the VHE range. Two flavors of the microquasar model 
can be found in the literature depending on whether 
hadronic or leptonic jet matter dominates the emission 
at such an energy range. Among leptonic jet mod- 
els, there are inverse Compton jet emission models in 
which X-rays and 7 -rays result from synchrotro n self- 
Compton processes ( Atovan fc Aharonian 19991 ). or in 
which the seed photons come from external sources, 
i . e companion star and/or accretion disk photons 



( Kaufman Bernado et al.l l2002l : iGeorganopoulos et al, 
2OO2I) . In the hadronic scenario, the gamma-ray emis- 



sion arises from the decay of neutral pions created 
in the inelastic collisions between relativistic protons 
in the jet and either the ion s in the stellar wind of 
the massive companion star ( Romero et al.l [2003h or 
nearby high-density regio ns (i.e. molecular clouds) 
( Bosch- Ramon et al.ll2005r ). 

Alternatively, relativistic particles can be injected 
in the surrounding medium by the wind from a young 

•^The Be star HD 215227, likely counterpart of the gamma-ray 
source AGL J2241 -1-4 454, has been suggested as a new candidate 
jWiliams et al.ll201(]|) . 



pulsar. In the pulsar wind scenario, the rotation of a 
young pulsar provides stable energy to the nonthermal 
relativistic particles in the shocked pulsar wind mate- 
rial outflowing from the binary companion. As in the 
microquasar-jet models, the 7-ray emission can be pro- 
duced by inverse Compton scattering of the relativis- 
tic particles from the pulsar wind on stellar photons 
( Torresllioiol and references therein). In the pulsar 
wind scenario the resolved radio emission is not due 
to a relativistic jet akin to those of microquasars, but 
arises instead from shoc ked pulsar w ind material out- 
flowing from the binary (iDubuslbOod) . 



1.2.4 "f-Cas like objects 

A growing number of early Be stars discovered in X- 
ray surveys exhibit strong X/optical flux correlations 
and X-ray luminosities intermediate between those of 
normal stars and those of m ost BeXB in quiescence 
( Lopes de Oliveira et al. 20061) . The optical properties 
are very similar to those of BeXB: i) the spectral type 
is always in the range BO-Bl Ill-Ve, and ii) show Ha 
equivalent widths stronger than —20 A. However, they 
differ from the typical BeXB in their X-ray properties: 
i) they show harder X-ray spectra that can be best 
fitted by a thin thermal plasma with T ^ 10* K, rather 
than a power law as seen typically in HMXBs, ii) there 
is no evidence for coherent pulsations in any of these 
systems but strong variability on time scales as short 
as 100 s is usually observed, and Hi) they do not exhibit 
large X-ray outbursts. 

The prototype of this group of sources is 7 Cas. Two 
models have been put forward to explain this type of 
system: accretion onto a compact object (most likely 
a white dwarf) and magnetically heated material be- 
tween the photo sphere of the B star and the inner 
part of its disc (jRobinson et al.l 120021 ). In the mag- 
netic corona model the hard X-rays result from high- 
energy particles that are emitted due to magnetic recon- 
nection, while the optical variability is due to changes 
in the density structure of the inner disc as a conse- 
quence of turbulence generated by changes in the mag- 
netic field. In support of this scenario there is the 
fact that X-ray fluxes show random variations with or- 
bital phase, thereby contradicting the binary accretion 
model, which predicts a substantial modulation. In 
favour of the binary model is the similarity of the X-ray 
spectra with those of cataclysmic variables, the rapid 
variability and large numbers of Be + white dwarf sys- 
tems predicted by the theory of evolution of massive 
binaries. 

In the remaining part of this report I shall concen- 
trate on the properties and variability of BcXBOnly 
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Fig. 4 IR excess and Ha emission are tiie two main observational characteristics of Be stars. 



when for the sake of the discussion a comparison with 
the behaviour of other type of X-ray binaries may be 
illustrative, will these other type of binaries be men- 
tioned. 



2 Optical/IR properties of BeXB 

The optical and infrared flux of a BeXB is completely 
dominated by the Be star companion. While the X-ray 
emission from BeXB provides information on the physi- 
cal conditions in the vicinity of the compact object, the 
optical and infrared emission reveals the physical state 
of the mass donor component. Since the fuel that pow- 
ers the X-ray emission, namely matter from a powerful 
stellar wind or from a equatorially concentrated denser 
disc, comes from the massive companion, observations 
of BeXB in the optical/IR are crucial to understand 
the physical conditions under which the neutron star is 
accreting. 

The two main observational characteristics of Be 
stars are the emission spectral lines, as opposed to the 
normal absorption photospheric lines and an excess of 
IR emission. Both properties, line emission and IR ex- 
cess, originate from extended circumstellar envelopes of 
ionised gas surrounding the equator of the B star. They 
result from free-free and free-bound emission from the 
disc, i.e. recombination of the optical and UV radiation 
from the central star ( Woolf et al. IQTOt Gehrz et al. 
I974I ). The assumption of a common origin for the 



on the geometry (size, shape) and dynamics (yeloc- 



line emission and IR excess is strongly supported by 
the correlation between the intensities of continuous IR 
emission, as measured either by colour indices (J — M), 
(J — K) or fluxes at a certain wavelength and the in- 
tensity of the H a line, as measure d either as equivalent 



width or fluxes (|Dachs fc Wamstcker 19 82: D achs et al 



19881 : iNeto fc de Freitas PachecoTllQsA . By studying 
the line variability one can obtain important constraints 



ics ( • 
ichiki 



ity arid density la. ws) of the envelope (jHanuschi 
199,4 lDachslll99i . 



1986 



Figure U] illustrates graphically these two properties 
for the system X Persei. When an equatorial disc is 
present the near-IR emission exceeds that predicted 
by model atmosphere (shown is a Kuruzc model with 
Tcs = 30000 K and log g = 3.5) and emission lines are 
prominently seen. 

2.1 Ha line profiles 

The Ha line is the prime indicator of the circumstel- 
lar disc state. Ha emissio n lines can be morphologi- 
cally divided in t wo classes (IHummel fc Vrancken|[l995i 
Silai et aLllioToh : class 1 are symmetric and includes 



symmetric double-peak, wine-bottle and shell profiles 
and class 2, which are asymmetric and show variability 
on time scales of years. Each class does not refer to dif- 
ferent groups of sources because individual Be stars can 
change from one to the other. These changes are nor- 
mally slow (of the order of years to decades in isolated 
Be stars and months to years in BeXB). Symmetric 
profiles are bel ieved to be gen erated in quasi-Keplerian 



discs (see e.g. iHummell 11994) . while asymmetric pro 



Hanuschik 


I995I 


Hummel fc HanuschikI 


1997) 



Most BeXB show asymmetric split Ha profiles. The 
peaks adopt the names of the relative position of their 
central wavelengths. In a spectrum with monotonically 
increasing wavelength values, the left peak is known 
as the "blue" (or violet) peak, while the right peak is 
named the "red" peak. V/R variability refers then to 
the variation of the relative strength of the blue to the 
red peak. Therefore the V/R ratio, defined as the ratio 
of violet-side to red-side peak intensities above contin- 
uum in units of continuum intensity represents a mea- 
sure of the asymmetry of the line. A more convenient 
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quantity is the logarithmic of this ratio, \og{V/ R), be- 
cause in this case, positive values of log(y/i?) corre- 
spond to a blue-dominated profile and negative values 
to a red-dominated line. The top panel of Fig. [5] dis- 
plays the V/R v ariability of the BeXB LS I +61 235 
(|Reig et al.|[2000l ). 

V/R variability can be explained in terms of a non- 
axisymmetrical equatorial disc in which a one-armed 
perturbatio n (a zone in the disc with higher density) 
pa gates (Okazakill991 . 1996 . 1997 : Papaloizou et al. 



pro 



1992t ISavoniie fc Heemskerklll993l ). Double-peak sym- 
metric profiles are expected when the high-density part 
is behind or in front of the star, while asymmetric pro- 
files are seen when the high-density perturbation is on 
one side of the disc. More precisely, when the high- 
density part of the disc perturbation is located on the 
side of the disc where the rotational motion is directed 
away from us, we see enhanced red emission (V < 
R), while when the high-density part is moving toward 
the obs erver, blue-dominat ed profiles V > R are ex- 
pected (jTelting et al.lll994l) . For systems seen at high- 
inclination angle, the two symmetric cases can be read- 
ily distinguished since the central depression between 
the two symmetric peaks would be more pronounced 
(reaching or going beyond the continuum, the so-called 
shell profile) , when the perturbation is between the ob- 
server and the star. The reason is that shell profiles are 
thought to be due to a perspective effect, namely, when 
the line of sight toward the star prob es the equato 



rial d isc and self-absorption is produced (jRivinius et al 
20061) . If the density perturbation revolves around the 



star in the same direction as the material in the disc 
(prograde p recession), the i i the V /R sequence would be 
according to lTelting et al. ( 1994) : V = R (perturbation 
behind the star) — > V > R — V ^ R (perturbation 
in front of the star, shell profile) — !► V < R. 

One prediction of the model is that no changes in 
the slope of the infrared continuum are expected be- 
cause the V/R variations are not the result of changes 
in the radial gradient of the circumstellar gas. This is 
exactly the behaviour that it is found in LS I -t-61 235 
(bottom panel of Fig. [5]). While the individual infrared 
photometric bands changed (AJ sa AH w AK ^ 0.3 
magnitudes ) the infrared colo urs (e.g. J^K) remained 
unchanged ( Reig et al. 2000l ). 

V/R variability is also seen in other lines, like Hel 
6678A. Since helium lines are gen erated at smaller disc 
radii than the hydrogen line s ([Hummel fc Vrancken 
1995t IJaschek fc Jaschekl 120041 ) . the asymmetry of the 
Hel line profiles indicates that the internal changes of 
the disc are global, affecting its entire structure. 

V/R quasi-periods in BeXB range from 1-5 yr (Ta- 
ble [5]) and are shorter than those seen in isolated Be 
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Fig. 5 Top: V/R variability observ ed in the Ha line of 
LS I -1-61 235. From lReig et all (|2000l '). Bottom: Evolution 
of the V/R ratio, J magnitude, J — K colour and t he H q 
equivalent width in LS I -|-61 235. From lReig et all (|2000h . 



stars, which are fo und in the ran ge 2-11 years with an 
average of 7 years (jOkazakilll997() . 



2.2 EW(Hq!) vs infrared colours 

If the infrared excess observed in Be stars is due to 
the same processes as those responsible for Balmer 
line emission, namely absorption and subsequent re- 
emission of the optical and UV light from the underly- 
ing star in the circumstellar envelope, then a correlation 
between IR colours and the strength of the hydrogen 



ported for isolated Be stars 


( 


Dachs & Wamsteker 1982 


Neto & de Freitas Pacheco 




19821: Dachs et al. 


1988) 


and BeXB (Coe et al.lll994 


, 20051 although mixed with 



isolated Be stars) . Figure |6] shows the first such dia- 
gram made from BeXB only. It displays the IR colour 
index { J—K)o as a function of the Ha equivalent width, 
EW(Hq;). for the sources listed in Table[3l The infrared 
colours were corrected for interstellar extinction. Only 
contemporaneous data (when the time difference be- 
tween the IR and optical observations was less than 
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-0.4 



-0.2 



(J-K)„ 



Fig. 6 Correlation between the equivalent width of the Ha 
line and the infrared colour J — K. Different colours repre- 
sent different systems. This correlation implies a common 
origin for the infrared excess and the line emission, namely, 
the circumstellar disc round the Be star. 



one month) were used. The second and third columns 
in Table |3] represent the amplitude of change in the 
infrared colour (J — K) and the EW(Ha) over the ob- 
served range, whereas Fig. [5] plots the actual values. 

In Fig. [S] we have included only those sources for 
which we are confident that the colour excess E{B — V) 
is free of circumstellar effects. Due to the surround- 
ing envelope, the use of photometric magnitudes and 
colours to derive astrophysical parameters, like E[B — 
V), may be misleading because they are contaminated 
by disc emission. The effect of the disc is to make 
the photometric indices to appear redder than a non- 
emitting B star of the same spectral type. The disc 
emission makes a maximum contribution to the op- 
tical (B — V) colour of a few tenths of a magnitude 
i Dachs et ailllQSsHHowells et al.ll200l[ ). The values of 
E{B — V) used to produce Fig. [6] were obtained from 



Table 3 List of BeXB used in Fig. [6] with colour excess 
and amplitude of change in (J — K)o and EW(HQf). 



X-ray source 


E{B - V) A{J-K)o 


AEW{Ha) 




(mag) 


(mag) 


(A) 


4U 0115+63 


1.65 


0.8 


20 


RX J0146. 9+6121 


0.93 


0.7 


10 


V 0332+53 


1.88 


0.4 


6 


4U 0352+30 


0.39 


1.1 


18 


RX J0440. 1+4431 


0.65 


0.3 


5 


lA 0535+262 


0.75 


1.0 


15 


4U 0728-25 


0.73 


0.1 


2 


RX J0812.4-3114 


0.65 


0.2 


5 


GRO J1008-57 


1.90 


1.0 


15 


lA 1118-616 


0.90 


1.5 


25 


4U 1145-619 


0.29 


2.8 


40 


EXO 2030+375 


3.8 


1.1 


15 


SAX J2103. 5+4545 


1.35 


0.1 


1 



either diffuse interstellar bands, disc-loss episodes or by 
disentar igling the circumstellar and interstellar redden- 
ing (see iFabregat fc Torrei6nlll998r ). 

Under these conditions, namely, simultaneity of the 
optical and infrared data and used of the extinction- 
corrected (J - K) index, the EW{Ha) - (J - K)o 
becomes a useful tool to estimate the extra reddening 
caused by the circumstellar disc. Once a value of the 
EW(Hq;) is known, one can look for the corresponding 
intrinsic (J — K) in Fig. [6l By comparing this value 
with that expected according to the spectral type, an 
estimated of the disc contribution on the index (J— K) 
can be obtained. The relatively large scatter, however, 
may limit the usefulness of the diagram. 

It should be noticed that although the Ha and in- 
frared emission correlate as expected for a common ori- 
gin in the disc, the spatial extension and precise lo- 
calisation in the disc are different. From lo ng-baseline 
interferometric observations in the K' band Giesetal 



(|2007i) found that the angular size of the infrared emis- 
sion is consistently smaller than that determined for 
the Ha emission. In other words, the near-IR emis- 
sion forms closer to the star than does the Ha emis- 
sion. Furthermore there is good evidence in s everal 



cases ( Clark et aLluOOlt iGrundstrom et al.ll2007[ ) that 



an increase in disc brightness occurs first in the near-IR 
flux excess and later in Ha as expected for an outwardly 
progressing density enhancement. This difference can 
be explained if the dependence with disc radius of the 
Ha optical de pth is less steep t han that of the infrared 
optical depth (jCies et al.ll2007t ). 



2.3 Identification of the optical counterparts 

With the improved sensitivities of the currently opera- 
tional space missions many new X-ray sources are being 
discovered. About one third of th e the X/7-ray sou rces 
in the 4th INTEGRAL catalogue |Bird et al.lboiol) are 
thought to be X-ray binaries, of which half are believed 
to contain early-type companions. An optical identi- 
fication is necessary to facilitate a complete study of 
these systems. Without a known counterpart, observa- 
tions are limited to X-ray energies, and hence our un- 
derstanding of the structure and dynamics of those sys- 
tems that remain optically unidentified is incomplete. 

If the only available information is that provided by 
an isolated X-ray detection then potential HMXBs are 
selected as those exhibiting properties of a magnetised 
neutron star, namely. X-ray pulsations and/or an ab- 
sorbed power-law continuum spectrum with an expo- 
nential cutoff at 10-30 keV and cyclotron absorption 
features. Extra information on the nature of the source 
can be obtained if the long-term X-ray variability is 
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B-V 



Fig. 7 Colour-colour diagrams and V-band images of tiie field around SAX J2103. 5+4545 (top) and IGR J01363+6610 
(bottom). The position of the confirmed optical counterpart in the colour-colour diagram is marked with a filled circle. 
Good candidates are those that occupy the upper part of the colour-colour diagram and lie inside or very close to the X-ray 
uncertainty region. For example, the stars that are close to SAX J2103.5-f 4545 in the colour-colour diagram lie far away 
from the satellite error circles. 



known (information provided by all-sky monitors, like 
ASM RXTE or BAT SWIFT). The presence of regular 
and periodic outbursts or unexpected giant outbursts 
reaching Eddington luminosity at their peaks may indi- 
cate the presence of a Be star. Likewise, erratic flaring, 
i.e., X-ray variability with changes in the X-ray inten- 
sity by a factor 3-5 over a few minutes might indicate 
the presence of a strong stellar wind from an early-type, 
probably evolved, companion. 

However, a system with X-ray data only will remain 
in the category of suspected HMXBs until a confirmed 
optical identification is performed. The most obvious 
observational features to look for is Ha emission and 
near-infrared excess emission. While the detection of 
these two observational features does not guarantee 
that the source is indeed an OBe type star, it defi- 
nitely narrows the number of candidates to a handful 
of sources. 

The size of the uncertainty in the location of the 
X-ray source, the so-called error circle, determines the 
type of observational technique to use. If the error circle 
is small (a few arcseconds) then the number of visible 
stars in the region is expected to be small and it is pos- 
sible to perform narrow-slit spectroscopic observations 
and look for early-type. Ha emitting line stars. If the 
error radius is large (>1 arcminute), then it is likely to 
include a large number of sources and hence narrow-slit 
spectroscopy becomes impractical. In this case colour- 
colour diagrams can be used to identify good candi- 



dates. Furthermore, uncertainty regions are given to 
a certain confidence level, hence it is possible that the 
true optical companion is located close but outside the 
X-ray error circle. 

Photometric detection of Be stars can be performed 
by selecting colours directly related to the excess emis- 
sion both in the Ha line (" red" colour) and the contin- 
uum ("blue" colour). Observing through a narrow filter 
centred on the Ha line and a wider filter also contain- 



ing this line, e.g. Johnson R filter (jReig et al.l 120051 ) 



Sloan-r ( Gutierrez- Soto 20061) would account for the 
emission line, while Johnson B — V (IReig et al.l 120051 ) 
or Stromgren b — y (jGutierrez- Sotoll2006l) can be used 



as reddening indicators. The use of Stromgren b, y, and 
narrow-band Ha photometry through (b — y, y) colour- 
magnitude and {b ~ y ~ Ha) colour-colour diagrams 
is pa rticularly suitable to identify Be stars in op en clus- 
ters (jCrebel et al.lll992l : iMcSwain fc Giesll2005[ ). 

Stars with a moderately or large Ha excess can be 
distinguished from the rest because they deviate from 
the general trend and occupy the upper left parts of the 
diagram. Be star are expected to show bluer colours, 
i.e., low {B — V) or {b — y), because they are early-type 
stars (although they normally appear redder than non- 
emitting B stars due to the circumstellar disc) and also 
larger (i.e., less negative) R — Ha colours because they 
show Ha in emission. Once the number of candidates is 
reduced to a few systems, narrow-slit spectroscopy be- 
comes feasible. The final step is to obtain a spectrum 



11 




Fig. 8 Top left: Q ~ K diagram for three well-known systems. BeXB systems tend to occupy the left-most part of this 
diagram. Top left: {b — y,y) diag ram of the open c l uster NGC 3766. Be stars (diamonds) are distinguished from all other 
stars (dots) in the diagram. From IMcS wain fc GiesI |200^. Bottom: IR photometric diagrams of the field around the BeXB 
lA 1118-615 (red cross). BeXB candidates (open circles) show up below the sequence of absorption-line stars (black dots). 



m the 4000-5000 A, and perform spectral classification. 
Figure [7] shows the V-band images and the correspond- 
ing {B — V, R — Ha) colour-colour diagrams of two 
recently identified systems 

At other wav elengths, similar photometric d iagrams 
have been used. iNegueruela fc Schurch ] (l2007t) defined 
the reddening-free quantity Q = {J—H) — 1.70{H—Ks), 
and created Q /Kg diagrams to separate early type from 
late type stars. The majority of stars in Galactic fields 
concentrate around Q=0.4-0.5, corresponding to field 
K and M stars, while early-type stars typically have 
Q ~ 0. The top left panel in Fig. [5] shows one such 
diagram with three well-kn o wn Be XB. 

More recently, iNesDolil (I2OIOI) built up infrared 
colour diagrams {Brj — Ks) — {H — K^) and {Hel — 
Ks) — [H — Kg) to identify emission-line objects. Br7 is 
the most prominent feature in Be stars in the K band, 
while He I 20580 A is found in early type Be stars, 
up to B2.5 IClark fc Steele! I2OOOI) . Emission-line ob- 
jects show up below the absorption-line stars sequence. 
Since Be stars present moderate emission line strengths, 
when compard to other groups (young stellar objects, 
cataclysmic variables, planetary nebulae). Be star can- 
didates arc those below the main sequence but that do 



not differ more than ^0.3 mag (Fig. The top right 
diagram in Fig. [5] also shows an example of the use of 
narrow-band filters to i dentify Be stars in clusters (from 
McSwain fc Giesll2005l) . 



3 X-ray properties 

Tradiationally, BeXB have been considered as transient 
X-ray binaries harbouring a fast rotating neutron star 
that orbits a Be stype star in a rather eccentric orbit 
(e > 0.3). Although X-Per was recognized as a BeXB 
since the beginning of the suggestion of the existence 



of th is type of objects in the early 70's (iMoffat et al 



19731) . it was the only BeXB showing permanent low 
level X-ray emission. As the number of BeXB sys- 
tems began to grow, the X-ray behaviour increased in 
complexity and new persistent sources were discovered 
iReig fc Rochelll999l ). 

Because all BeXB are variable and the detectability 
of a source depends on the sensitivity of the detector 
used, the distinction between transient and persistent 
sources undergoing a large intensity increase is some- 
times ambiguous. The term transient is normally used 
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Fig. 9 Long-term light curves of EXO 2030+375 (left) and 4U 0115+63 (right). Although the two systems have been 
seen to display both type I and type II outbursts, type II outbursts in EXO 2030+375 are rare whereas type I are almost 
always present. In 4U 0115+63 the opposite occurs. 



when the variability of the source exceeds, at least, two 
orders of magnitude with respect to the quiescent state, 
which may correspond to a non-detection state. 

3.1 Persistent versus transient X-ray emission 



Reig fc Roch3 |l999l ) proposed the following character- 
istics differentiating between persistent and transient 
BeXB. Persistent sources are characterised by 

• low X-ray luminosity, L2_20koV ^ 10^4-35 gj-g g-i_ 

• relatively quiet systems showing flat light curves with 
sporadic and unpredicted increases in intensity by 
less than an order of magnitude. 

• slowly rotating pulsars, Pgpin ^ 200 seconds. 

• absence of, or very weak, iron line at ^6.4 keV, in- 
dicative of only small amounts of material in the 
vicinity of the neutron star. 

These X-ray properties could be accommodated 
in systems with wide (Porb ?J 200 days) and low- 
eccentric orbits (e < 0.2). A thermal excess of black- 
body type, with high temperature (kT > 1 keV) and 
small emission area {R < 0.5 km), has recently been 
suggested as another common featur e of this type 
of sources ( La Palombara et al. I boool and references 
therein). Members of this group are: X Per, RX 
J0146.9+6121/LS I +61 235, RX J0440.9+4431/BSD 
24-491 and RX J1037.5-564/LS 1698. 

In contrast, the X-ray behaviour of transient BeXB 
is characterised by two type of outbursting activity: 

• Type I outbursts. These are regular and periodic 
(or quasiperiodic) outbursts, normally peaking at or 
close to periastron passage of the neutron star. They 
are short-lived, i.e., tend to cover a relatively small 
fraction of the orbital period (typically 0.2-0.3 Porb)- 
The X-ray flux increases by about one order of magni- 
tude with respect to the pre-outburst state, reaching 
peak luminosities < 10^^ erg s~^. 



• Type II outbursts represent major increases of the 
X-ray flux, 10^ — 10^ times that at quiescence. They 
reach the Eddington luminosity for a neutron star 
and become the brightest objects of the X-ray sky. 
They do not show any preferred orbital phase and 
last for a large fraction of an orbital period or even 
for several orbital periods. The formati on of an ac 



cretio i i disc during Type II outbursts (IKriss et al 



1983 



Motch etal 



19911: iHavasaki fc Ozakakil 12004 : 
Wilson et al.l |2008[ ) may occur. The discovery of 
quasi-periodic oscillations in some systems (?, and 
reference therein) would support this scenario. The 
presence of an accretion disc also helps explain the 
large and steady spin-up rates seen during the giant 
outbursts, which are difficult to account for by means 
of direct accretion. 

Figure [5] shows these two types of X-ray variabil- 
ity. The long-term X-ray emission of EXO 2030+375 
appears clearly modulated by the orbital period of 46 
days just before and after the type II outburst at MJD 
53950. 4U 0115+63 has exhibited four type II outburst 
in the time interval 1996-2009. 

Type II outbursts are major events. Since the fuel 
that powers these outbursts comes from the material in 
the circumstellar disc, one would expect major changes 
in the structure of the disc. The disruption should be 
observable in the parameters that best sample the phys- 
ical condition in the disc, namely the Ha line param- 
eters (strength and shape) and the IR colours and in- 
dices. Figure [TUl displays the long-term evolution of the 
optical/IR and Ha equivalent width of 4U 0115+63. 
Vertical dotted lines mark the occurrence of type II X- 
ray outbursts. As it can be seen, the large amplitude 
changes in the optical and IR photometric and spec- 
troscopic parameters after the X-ray outbursts indicate 
major disruption in the physical conditions of the cir- 
cumstellar. In fact, the type II outbursts led to the 
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Fig. 10 Long-term evolution of the optical/IR magnitudes and the Ha equivalent width and pea k separation of d ouble- 
peak profiles of 4U 0115H-63 . Dotted lines mark the occurrence of X-ray outbursts. Adapted from lReig et alj (|2007l ') 



complete disappearance of the disc. Evidence for the 
loss of the disc stems from the Ha equivalent width, 
which changed sign from negative, indicating that the 
spectral line was in emission before the outburst to pos- 
itive, that is absorption, after the outburst. 

However, type II outbursts are not always followed 
by a disc-loss phase. Figure [11] shows the case of 
lA 0535-^262 and KS 1947-^300. In the case of lA 
0535+262, the Ha equivalent width decreased signif- 
icantly after the X-ray outburst, indicating a major 
disruption in the structure of the disc, but the Ha 
equivalent width remained negative. The case of KS 
1947+300 is hard to explain. Apparently, the X-ray 
outbursts produced no effect on the disc. The Ha equiv- 
alent width did not change significantly after the out- 
burst. 

3.2 Periodic variability: X-ray pulsations 

Observations of X-ray binaries show considerable vari- 
ability on a wide range of time scales in all wavelengths. 
In LMXB, in the so-called millisecond pulsars, the X- 
ray periodic variability goes down to less than a mil- 
lisecond. In BeXB the shortest timescales are of the 
order of few second^, and correspond to the X-ray pul- 
sations. 



^There are at least two systems, A0535-66 and SAX J0635+0533, 
with B-type companions that exhibit millisecond pulsations. 
However, they are thought to be rotation-powered pulsars. 



The detection of coherent pulsations from an accret- 
ing X-ray source provides one of the strongest evidence 
that the compact object is a neutron star. The pulses of 
high-energy radiation are due to a misalignment of the 
neutron star's rotation axis and its magnetic axis. Gas 
is accreted from the stellar companion and is channeled 
by the neutron star's magnetic field onto the magnetic 
poles producing two or more localized X-ray hot spots. 
These hot spots move into and out of view as the neu- 
tron star spins, giving rise to regular X-ray pulses. 

Almost all BeXB show X-ray pulsations. The pul- 
sation period of an X-ray pulsar with a Be compan- 
ion can be as short few seconds (4U 0115+63) 
and as long as few tens of minutes (LS I +61 235). 
The longest periods in main-seq uence HMXBs are those 
of 4U 2206+54 with 1.54 hr faeig et all l2009h and 



200' 



How- 



IGR J01583+6713 with 5.47 hr (jWand 
ever, these two sources are somehow peculiar and can- 
not be considered as canonical BcXB. Although 4U 
2206+54 exhibits a double-peaked emission Ha profile 
and the optical companion is an 09.5 main-sequence 
star, its optical spectrun i is complex and far from that 



of a "typical" Be star |Blav et all l2006l ). The 5.47- 



hr period of IGR J01583+6713 needs confirmation as 
a shorter (469 s) periodicity has also been suggested 
(jKaur et al.ll2008( ). Among HMXBs with a supergiant 
companion, 2 S 0114+65 shows t he longest spin period 
with 2.78 hr |Finlev et al.lll994l) . The only BeXB for 
which no pulsations have been detected is LS I +61 



14 



2000 2001 2002 2003 2004 2005 2006 2007 




52000 52500 53000 53500 54000 54500 
MJD 




,< -18- 



2000 2001 2002 2003 2004 2005 2006 2007 2008 

I ' I ' I I I ' I r 

KS 1947+300 



I , I , I , I , I , I 

52000 52500 53000 53500 54000 54500 
MJD 



Fig. 11 Long-term evolution of the X-ray intensity (1.3-12.1 keV) and the Ha equivalent of lA 0535+262 (left) and KS 
1947+300 (right). 



303/2E 0236.6+6101. The nature of this system is un- 
der debate. The lack of pulsations can be explained if 
it contains a black hole, that is, if LS I +61 303 is a 
microquasar. Radial vel ocity studies do not co mpletely 
rule out this possibility (jCasares et al.ll2005al) . On the 
other hand, radio observations favour the presence of a 
neutron star. However, in this case LS I +61 303 would 
contain a rotation-powered pu lsar with nebula and not 
an accretion-powered pulsar (jPhawan et al.l l2006l . see 
also Sect. [1X31 . 

During major outbursts, transient accretion-powered 
pulsars show spin-up episodes, with the spin-up rate 
increasing with accretion luminosity (Fig. [T^ . This is 
consistent with the simple model o f disc magnetosphere 
interaction (jChosh fc Lamblll979l ). in which, as the ac- 
cretion rate increases, the rate of angular momentum 
transfer from the accretion flow to the neutron star in- 
creases. At the limit of very low mass accretion rate, 
the neutron stars are expected to spin down, because 
the magnetic field lines that couple to the outer, slower 
accretion flow remove spin angular momentum from the 
neutron star. 

Spin periods in the range 1-1000 s can be ex- 
plained according to the current understanding of 
the neut ron star's sp i n evo l ution in a close binary 



called the subsonic p ro peller phase (jPavies fc Pringle 
19811 : IlkhsanovlliooTal lbl. l2007h . 



system IIDavies et al 



1979 



I Davies 
20061) 



Pringle 1981 



Zhang et all 12004 iDai et al.l 120061 ). Longer periods 
are, however, difficult to accommodate in the standard 
theory. The spin evolution of a neutron star in a bi- 
nary system is divided into three basic phases or states. 
Each state is characterised by a different energy release 
mechanism and corresponds to a different evolution- 
ary stage of the neutron star. These phases are known 
as the pulsar (or ejector) phase, the propeller phase 
and the accretion phase. To explain the existence of 
long-period pulsars (Pspin ^ 1000 s), a subphasc of the 
propeller state needs to be invoked. This subphasc is 



After a supernova explosion a neutron star is formed 
as a rapidly rotating radio pulsar. The spin period at 
birth is of the order of a fraction of a second and its 
magnetic field of the order of lO^'^ G. In this phase, the 
neutron star radiates at the expense of its rotational 
energy. 

The condition for the transition of a slowing-down 
neutron star from the ejector phase into the supersonic 
propeller phase occurs when the pressure of the rel- 
ativistic wind ejected by the neutron star (magnetic 
pressure) is balanced by the pressure of the plasma sur- 
rounding the neutron star (ram pressure) at the grav- 
itational capture or accretion radius race = '^GAIx/v'^, 
where is the neutron star mass and is the stel- 
lar wind velocity. The plasma will tend to accrete onto 
the neutron star under the action of gravitation. The 
electromagnetic field, however, will obstruct this pro- 
cess, and the accreting matter will come to a stop at 
the magnetospheric boundary. The accreted material 
cannot penetrate any further because the drag exerted 
by the magnetic field is super-Keplerian. The in-falling 
matter is accelerated outward, taking away the angular 
momentum of the neutron star. Thus the spin-down 
of the neutron star continues until the centrifugal and 
gravity forces acting on a particle that corotates with 
the neutron star balance. This occurs at the so-called 
equilibrium period. 

If the neutron star spin period is larger than the equi- 
librium period, then the centrifugal barrier ceases to be 
effective and plasma that penetrates into the magneto- 
sphere is able to reach the neutron star surface, moving 
along the magnetic field lines in the direction of the 
magnetic poles. 

Since the material accreted carries substantial angu- 
lar momentum, the neutron star will experience strong 
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spin-up. Therefore, in principle, the equilibrium period 
given by 



Mr, 



-5/7 



represents the maximum period for a given mass accre- 
tion rate. /X30 and M15 are the magnetic moment and 
mass accretion rate in units of 10^° G cm^ and 10^^ g 
s~^, respectively. 

Pcq ^ 10"^ s can be achieved if the magnetic field 
strength in units of 10^^ G, B12 > 100 or if M15 < 
10~^. None of these conditions is supported by the 
observations. The condition of supercritical magnetic 
fields would imply that the neutron star is a magnetar. 
However, the detection of cyclotron lines in the range 
10-100 kcV indicates that the magnetic field strength 
is of the order a few times 10^^ q (^^^ - 1 - 6). The 
low mass accretion rate condition would imply an X-ray 
luminosity two to three orders of magnitude lower than 
observed. 

In order to explain spin periods longer than ~ 10^ 
s, another phase in the spin-period evolution has to 
be considered. This phase is called subsonic propeller 
phase and occurs p rior t o the accretion-powered phase 
(|lkhsanov|[2001al lbl. F 



20071) . 



We indicated above that if Pspin > ^cq then matter 
can be accreted onto the neutron st ar surface. How- 
ever, t his is actually not always true. Davies fc Pringlj 
(jl98lh showed that the star's magnetosphere dur- 
ing the propeller spin-down epoch is surrounded by 
a hot, spherical quasi-static envelope. This enve- 
lope extends from the magnetospheric radius r™ = 
/ . \ 2/7 

( /M \J2GMx I up to the accretion radius race = 



2GMx/v'^. The interaction between the magnetosphere 
and the envelope leads to the deceleration of the rota- 
tion rate of the neutron star: the rotational energy loss 
by the star is convected up through the envelope by the 
turbulent motions and lost through its outer boundary. 
The neutron star remains in the propeller state as long 
as the energy input to the envelope due to the propeller 
action by the star dominates the radiative losses from 
the envelope plasma. 

Steady accretion under the condition Pgpin > Pcq can 
be realized only if the cooling of the envelope plasma 
(due to radiation and convective motions) dominates 
the energy input. The subsonic propeller phase corre- 
sponds to the phase when this cooling occurs. 

The condition for accretion then translates to Pspin > 
Pbr, wher e Pbr is the so -called break period, which is 
given by (|lkhsanovll2007l) 



Pbr 
Pbr 
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While Pspin < Pbr the temperature of the envelope re- 
mains higher than the free-fall temperature (Tff(r„j) = 
GMj-mp/krm, where rup is the proton mass and k the 
Boltzmann constant) and no accretion is possible. 




Flu* (la-^ erg cm'^ 5"') 



GRO J1744-28 ^ 



Fig. 12 Observed relationships between flux and pul- 
sar spin-up rate Left: A 0535-1-262 durmg the 1994 gi- 
ant and GRO J1744 -I-28 during the 1995-1996 outburst 
(|Bildsten et"al]|l997l). Right EXO 20 30+375 during the 
2006 giant outburst (jWilson et al.ll2008l ) 



3.3 Aperiodic variability 

Aperiodic variability refers to the variability caused by 
proce sses that do not repeat t hemselves in a periodic 

20061) . Pulsations, eclipses. 



van der Klis 1995 



way 

outbursts are excluded from the analysis of aperiodic 
variability studies. In contrast, flickering, irregular flar- 
ing, fluctuations are the subject of study. The main 
tool used for studying the aperiodic variability is the 
Fourier power spectrum, or a Power Spectral Density 
(PSD), of the count-rate time series. A number of vari- 
ability components, or power-spectral components, to- 
gether make up the power spectrum. An aperiodic com- 
ponent by definition covers several frequency resolution 
elements. Broad features are called noise and narrow 
features are called quasi-periodic oscillations (QPOs). 
The so-called Q-factor determines whether a noise com- 
ponent deserves the name QPO. If the Q-factor, defined 
as the ratio between the characteristic frequency of the 
component over its width, Q = v/FWHM, is larger 
than 2, then the noise component is considered to be 
narrow enough to be called a QPO. 

Aperiodic variability studies are normally carried out 
in correlation with the position of the source in the 
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Fig. 13 Hardness-intensity diagrams traced by various BeXB during giant (type II) outbursts observed by RXTE. The 
hard colour, HC, is the ratio between the count rate in the energy range 15-30 keV over 10-15 keV. 



colour-colour diagram and/or hardness-intensity dia- 
gram. An X-ray color is a 'hardness' ratio between 
the photon counts in two broad bands. A colour-colour 
diagram (CD) is the plot of one color vs. another cal- 
culated in different energy bands. The higher energy 
range color is called hard colour, while the lower energy 
range is the soft colour. A hardness-intensity diagram 
(HID) is instead the plot of a colour vs. the count-rate 
in some broad spectral range. Figure [T51 shows the HID 
of eight BeXB during type 11 outbursts. The CD/HID 
are very useful tools to investigate the spectral varia- 
tions of an X-ray source because they are model- and 
instrument-independent and reflect the intrinsic prop- 
erties of the system. 

Together with the functional form of the variability 
components ("noise"), CD and HID have revolutionised 
our understanding of X-ray binaries, unveiling a whole 
new phenomenology on the spectral and timing prop- 
erties of X-ray binaries by introducing the notion of 
source states. A state is defined by the appearance of a 
spectral (i.e. power-law, blackbody) or variability com- 
ponent (i.e. Lorentzian) associated with a particular 
and well-defined position of the source in the CD/HID. 
A transition between states takes place when the rela- 
tive strength of the spectral or variability components 
varies an d the source motio n in the CD/HID changes 
direction ( van der Klij 2006 ). 



In LMXBs and BHBs the aperiodic variability is 
thought to originate in the irregular nature of the in- 
ner accretion fiow, as an accretion disc is assumed to 
exist in this type of systems. In HMXBs the details of 
how the magnetic field affects the accretion flow are not 
fully understood. 

While there are numerous references in the litera- 
ture on the application of CD/HID analysis on LMXBs 
and BHBs, very little work of this typ e has been done 
on HM XBs |Belloni fc Hasinger][l99o[) . Recently, |Rei3 
(|2Q08l) performed the flrst systematic analysis of the 
spectral and rapid aperiodic variability of Be/X-ray bi- 
naries during type II outbursts. He found similarities 
and differences with respect to the behaviour of LMXBs 
and BHBs. 

The similarities include: i) existence of spectral 
branches in the CD/HID. At high and intermediate 
flux the sources move along a diagonal branch (DB), 
while at very low flux the soft colour decreases while 
the hard colour remains constant defining a horizontal 
branch (HB); ii) the source does not jump through the 
diagram but moves smoothly, i.e. without jumps, fol- 
lowing the pattern; Hi) the broad noise components in 
the HB are more variable (in terms of fractional rms) 
than in the DB; iv) the power spectra of HMXBs can be 
fitted acceptably with a sum of a small number (2-4) 
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of Lorcntzian functions (Fig. [T4|) . All sources present 
peaked noise in certain spectral states. 

While the differences are: i) presence of hysteresis 
patterns in the CD/HID that might be a manifestation 
of the magnetic field (through cyclotron resonant scat- 
tering features); ii) slower motion along the spectral 
branches, hours to days in LMXBs, days to weeks in 
HMXBs; in) the characteristic time scales implied by 
the noise components are about one order of magnitude 
longer in HMXBs. There are about 14 H MXBs pulsars 



that show QPOs in their power spectra (IJames et al. 
The frequencies of these QPOs range from 1- 



200 mHz. In contrast, in LMXBs and BHBs QPO range 
1-1000 Hz; iv) the so-called L5 noise, whose central fre- 
quency coincides with the frequency of the fundamental 
peak of the pulse period and suggests a strong coupling 
between the periodic and aperiodic noise components. 
This coupling implies that the instabilities in the ac- 
cretion flow that give rise to the aperiodic variability 
must travel all the way down to the neutron sta r sur- 
face iLazzati k Stellal Il997t iBurderi etliD Il997t ). and 
v) no apparent correlation between the power spectral 
parameters (characteristic frequencies, rms) and mass 
accretion rate. 



4 Circumstellar disc 

4.1 Disc formation 

Be stars are fast rotators. They have, on average, larger 
observed rotatio nal velocities than B stars as a group 
(jSlettebaklll982i ). The determination of the rotational 
velocity is believed to be a crucial parameter in the 
formation of the circumstellar disc. A rotational veloc- 
ity close to the break-up or critical velocity (i.e. the 
velocity at which centrifugal forces balance Newtonian 
gravity) reduces the effective equatorial gravity to the 
extent that weak processes such as, gas pressure and/or 
non-radial pulsations may trigger the ejection of photo- 
spheric matter with sufficient energy and angular mo- 
mentum to make it spin up into a Keplerian disc. 

At present, there is no consensus on how close to 
critical velocity Be stars rotate, nor on how Be stars 
reached such high spinning rates. Observations suggest 
that a large fracti on of Be stars r otate at 70-80% of 
the c ritical value (ISlettebaklll98i Porter] Il996i 1^ 
200 ll ). However, iTownsend et al.l |2004l ) claimed that 



if gravity darkening is taken into account then most 
or all Be stars may be rotating very nearly at critical 
velocity. Whether the rotational velocities of Be stars 
reflect the initial distribution of angular momentum or 
they have been spun up as a result of binary evolution 



is still a matter of debate. There is growing evidence 
that supports the idea that some Be star s were spu n 
up by mass transfer in interacting binaries ( GiesI 2000t ) . 
The former mass donors are now found as neutron stars 
(i.e. BeXB) in some cases and as a small hot subdwarfs 
stripped of its outer envelope (Phi Per, 59 Cyg, FY 
CMa ) in other cases (IGies et al.l Il998t iRivinius et al 



l2004t iMaintz et al.lboost IPeters et al.ll2008l) 

The projected rotational velocities are determined 
by measuring the wid th of certain spectral lines (see 
e.g. Steele et al. 19991) . as it is generally assumed that 
the width of the spectral lines and the rotational ve- 
locity are linearly related. This is true even with the 
inclusion of limb and intrinsic darkening, but if gravi- 
tational darkening is taken into account, then for stars 
rotating above ^80% of critical speeds the linear rela- 
tionship is no longer applicable. Any increase in the 
rotational velocity is accompanied by alrnost no change 
the line width ( Collins fc TruaxI 1995 : Cohen et al 



20051) . As a result, the observational parameter vsini 
may systematically underes timate the true proje cted 
equatori al rotatioii veloci ty ( Townsend et al. 20041 ) . In 
contrast, ICranmeii (|2005() found that this may be true 
for the late- type Be stars (i.e., spectral types B3 and 
later), but early- type Be stars do seem to be consistent 
with a range of intrinsic rotation speeds between ^40% 
and 100% critical. For stars rotating sufficiently below 
the critical velocity, any physical model for the origin of 
Be-star discs requires a significant increase in angular 
momentum above the photosphere. 

While the origin of the gas in these discs is agreed 
to be material ejected from the stellar photosphere, the 
precise mechanism triggering this mass ejection remains 
elusive. One of the first models that tried to explain the 
formati on of circumstellar discs is the wind-compressed 
model (jBjorkman fc Cassinellil 119931 ). The fundamen- 
tal idea behind this model is that in a rotating wind, 
the material in the wind tends to orbit the star dur- 
ing the time it is accelerated outward. If the outward 
acceleration is small compared to the rotation, then 
the radiation-driven wind flows toward the equator. 
Complimentary streamlines from opposite hemispheres 
would intercept each other at the equatorial plane giv- 
ing rise to a pair of shocks above and below the equator 
and eventually the disc. Note that the disc is not a re- 
sult of enhanced equatorial mass-loss. A major problem 
of this model arises when non-radial components are 
included in the radiation force, as these components 
lead to an effective suppression of the equatorward flow 
needed to compress the wind IjOwocki et al.lll996l ). In 
addition, this mo del was unab le to reproduce the IR 
excess of Be stars (jPorterlll997l ). 



The magnetically torqued disc model (iBrown et al 



19981) extends the ideas of the wind compressed disc 
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Fig. 14 Power spectra of 4U 0115+63 at various stages o f the 2 004 g iant outburst. The number at the top left of each 
panel is the X-ray luminosity in units of 10^'^ erg s^^. From lReiS (|20Q8l ). 



model by adding magnetic steering and torquing up of 
the angular momentum of the wind. The up-flowing 
wind from near the star is magnet ically channeled back 
to the equator. This obviates the lOwocki et al. I (Il996l) 
objections such as the driving of a polar flow by the non- 
radial radiation field. The magnetically torqued disc 
model requires a dipole-like magnetic field to channel 
the flow of wind material into a disc region along the 
equatorial plane. 

An alternative scenario that has been put forward 
to elevate disc matter to orbits well above the stel- 
lar surface and which is increasing in popularity is 
non-radial pulsations (jCranmeii l2009l and references 
therein). Non-radial pulsations may give rise to out- 
wardly propagating circumstcUar waves, which inject 
enough angular momentum into the upper atmosphere 
to spin up a Keplcrian disc. The main uncertainty in 
this model is the connection between pulsations and 
resonant waves since pulsations themselves are evanes- 
cent in the stellar photosphere. However, once the 
resonant waves have been formed, they grow in am- 
plitude with increasing height, begin to propagate up- 
wards, and steepen into shocks. The resulting dissipa- 
tion would create substantial wave pressure that both 
increases the atmospheric scale height and transports 
angular momentum upwards. 

Once the disc is formed, the viscous dccretion disc 



model (|Lee et al.l llQQlt lOkazakil 120011) stands up as a 
very promising model that explains many of the ob- 
served features in discs around Be stars. In this model. 



matter supplied from the equatorial surface of the star 
gradually drifts outward by the viscous torque and 
forms the disc. Decretion discs operate in a simi- 
lar fashion to accretion discs, except that the sign of 
M (mass decretion/accretion rate) is opposite. Vis- 
cosity is treated following the Shakura-Sunyaev's a- 
prescription. The outflow is highly subsonic in the inner 
part of the disc and close to Keplcrian. Although as- 
pects of the modeling (e.g., the supply of matter with 
sufficient angular momentum to the inner edge of the 
disc or the interaction of the stellar radiation with the 
disc matter) are not fully solved yet, the model is able 
to explain the very low outflow velocities (the observed 
upper limits are, at most, a few km s^^) and also the 
V/R variability observed in Be stars. 

4.2 Disc-loss episodes 

The use of photometric calibration for the determina- 
tion of the spectral type and luminosity class in a Be 
star is not as straightforward as in a non-emission-line 
B-type star owing to the presence of the surround- 
ing envelope, which distorts the characteristic photo- 
spheric spectrum. Be stars appear redder than the 
non-emission B stars, owing to the additional redden- 
ing caused by the hydrogen free-bound and free-free 
processes in the circumstellar envelope. 

Thus, in a Be star, one has to correct for both cir- 
cumstellar and interstellar reddening before any cali- 
bration can be used. There is no easy way to decouple 
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Table 4 Rotational velocity of Be/X-ray binaries 





Opticcil 


Spccticil 


Inclination 


V sin z 




source 


counterpart 


type 


angle (°) 


(km s^-'-) 




4U 0115+634 


V635 Cas 


B0.2V 


43 


300+50 


1 


RX J0146.9+6121 


LS I +61 235 


Bim-v 




200+30 


2 


V 0332+53 


BQ Cam 


08-9V 


<10 


<150 


3 


X-Per 


HD 24534 


09.5III 


23-30 


215+10 


4,5 


RX J0440.9+4431 


LS V +44 17 


Bim-v 




235+15 


6 


lA 0535+262 


HD 245770 


09.7III 


28-35 


225+10 


7,8 


RX J0812.4-3114 


LS 992 


B0.5III 




240+20 


9 


lA 1118-615 


Hen 3-640 


09.5IV 




300 


10 


4U 1145-619 


V801 Cen 


B0.2III 


<45 


250+30, 290 


11,12 


4U 1258-61 


V850 Cen 


B2V 


90 sh 




13 


SAX J2103.5+4545 




BOV 




240+20 


14 



Dcleado-Martf et al. 


J200l|), 


Janot-Pacheco et al. 


Jl98ll). f 



booth. [21 iReig et al. lll997^ 



[61 iRcig ct al. l bOOSi) [7 1 



[31 iNegueruela et al.l lll999ll [ 41 iLyubimkov et al.l 
Haigh et all (200^) [8] iGrundstrom et al " 



teooih [10] 

19801). [141 iReig et al] i2004l) 



1997), [9] iReig et al 



Table 5 Disc- loss episodes in Be/X-ray binaries. 



Source 


Spectral 


Vo 


{B - V)„ 


EW(Ha) 


^disc 


^loss 


Reference 


name 


type 






(A) 


(year) 


(day) 




4U 0115+63 


B0.2Ve 


15.5 


1.43 


+2.4 


3-5 


800 


1,2 


X Per 


BOVe 


6.80 


0.14 


+2.7 


7 


800 


3,4 


lA 0535+262 


09.7IIIe 


9.44 


0.41 


+2.5 


4-5 


600 


5,6 


RX J0812.4-3114 


B0.2IVe 


12.5 


0.4 


+1.4 


4 




7 


RX J0440.9+4431 


B0.2Ve 






-0.5 


>10 


<200 


8 


ICR J06074+2205 


B0.5Ve 






+2.2 


>5 




9 


SAX J2103.5+4545 


BOVe 


14.35 


1.07 


+2.8 


1.5-2 


1000 


10 



Note: [11 INegueruela et al.l bOO^ij. [2] iReig et a l.' 1*2007*), [3] iRoche et aO lll997^. [ 41 IClark et al.] ll200ll) [51 iHaigh et all l|l999l ). [6] 
IClark et al.l lll998l) , [7] lReig et al.l bOOll) . [81 iRelg ct al.. (i200g) [91 iReig et al.l (120103) . [10] iReig et al.l bOlOal) 
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these two reddening contributions. Although some iter- 
ative procedure has been proposed to correct for both 
circumstellar and interstellar reddeni ng from narrow- 
band photometry ( Fabregat fc Torrei6n,1998 | ) or using 



infrared colours ( Howells et alj 120011 : iDoughertv et al 



19941 ) ■ the most reliable measurement of the true pho- 
tometric colours of a Be star is when the disc contri- 
bution from the disc is negligible. That is the reason 
that disc-loss episodes are so important. The main ac- 
cepted disc- loss indicator is the Ha line (Fig.|4]). When 
the star loses the disc, the Ha line shows an absorp- 
tion profile and the X-ray activity ceases or is largely 
reduced. 

Table [5] list the BeXB that have gone through disc- 
loss phases. It also shows the typical time scales as- 
sociated with disc variability and the intrinsic colours 
of the underlying B star (not corrected for interstel- 
lar reddening). Tdisc is the typical duration of forma- 
tion/dissipation of the circumstellar disc and Tioss rep- 
resents the maximum duration of the disc-loss phase. 



5 Observational evidence of the disc-neutron 
star interaction 

Be stars may also exist as single objects, i.e., not form- 
ing part of an X-ray binary system. There is general 
consensus that the optical/IR long-term (months to 
years) variability seen in Be stars is caused by the equa- 
torial disc around the massive star. In BeXB, the X- 
ray variability is attributed to changes in the accretion 
rate. Since the circumstellar disc constitutes the fuel 
that powers the X-ray machinery, the X-ray variability 
can also be ascribed to changes in the disc. The Be star 
in a BeXB is assumed to have the same physical prop- 
erties (mass, radius, luminosity) as isolated Be stars. 
Therefore, any difference in the variability patterns of 
isolated and BeXB must be the result of different struc- 
ture and physical properties of the disc (size, density) . 

The following questions arise: Are the equatorial 
discs of isolated Be stars different from those of Be stars 
forming part of an X-ray binary? If yes, what causes 
this difference? 

The obvious difference between an isolated Be star 
and a BeXB is the presence of a neutron star in the 
latter. Since the underlying Be star has the same phys- 
ical characteristics in both type of systems, we conclude 
that any difference in the properties of the disc in bi- 
naries must be attributed to the neutron star. 

Traditionally, the interaction between the two com- 
ponents in a BeXB was thought to go into one direc- 
tion, namely, from the massive companion to the neu- 
tron star, through the transfer of mass. Matter from 



the Be star captured by the neutron star is accreted 
and powers the X-ray emission. Given the large mass 
ratio of the two components it was believed that the 
neutron star exerted little effect on the massive com- 
panion. This picture is now changing, as new observa- 
tional evidence for the interaction on the neutron star 
and the circumstellar disc emerges. 

In this section, I summarise the observational ev- 
idence gathered in recent years of how the compact 
object affects the optical/IR emission properties of a 
BeXB. The line of reasoning is the following: BeXB and 
single Be show different optical/IR patterns of variabil- 
ity because the equatorial disc has different physical 
properties. The disc exhibits different properties be- 
cause of the presence or absence of an external agent, 
that is, the neutron star. The neutron star controls 
somehow the evolution of the disc because it prevents 
its free expansion and causes truncation. 

In searching for observational evidence for the effects 
that the neutron star produces on the circumstellar disc 
of the Be star in BeXB, I have followed two different 
approaches 

1. Congelations with the orbital period. Since BeXB are 
in general moderately eccentric systems, any prop- 
erty that correlates with the orbital period of the 
system can be interpreted as due to the action of 
the neutron star. The logic behind this idea is that 
the shorter the orbital period, the smaller the or- 
bit. Hence circumstellar discs in narrow-obit sys- 
tems must be more affected by the gravitational pull 
exerted by the neutron star during each periastron 
passage than systems with longer orbital periods. 

2. Comparison of the variability patterns in isolated Be 
stars and BeXB. If there is no physical difference be- 
tween the B star in an isolated Be and the B star in 
a BeXB, and the optical/IR variability is driven by 
changes in the disc, then any difference in the ob- 
servational properties of these two types of systems 
must come from the effect of the neutron star on the 
disc. The objective here is to search for such differ- 
ences in the evolutionary patterns of various physi- 
cal parameters (disc density and radius, the strength 
and shape of spectral lines, particularly Ha, and/or 
characteristic time scales). 

This section begins with a summary of the viscous 
decrction disc model (see also Sect. 14. ip and it is fol- 
lowed by a compilation of the observational evidence of 
the interaction of the neutron star with the circumstel- 
lar envelope of the Be star. 

5.1 Disc truncation 

The idea of disc trun c ation in BeXB was first sug- 
gested by iReig et al.l (|l997() from an observational 
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Fig. 15 Porb-EW(HQ) diagram. Milky way objects have 
been represented by circles and BeXB in the Small Mag- 
ellanic Clouds by squares. The linear regression for the 
combined data sample is shown with a solid line, while the 
dashed line represents a power-law fit with exponent n=4/3. 



point of view and the theory subsequently developed 
by Okazaki and collaborator s in a series of papers: 
Negueruela fc Okazakil (120011) : lOkazaki fc Negueruela 
(|200ll ): IOkazaki et all (|2002l) . 

The central idea behind the viscous decretion disc 
models is that Be stars' discs are supported by viscos- 
ity. In this respect, the viscous excretion discs are very 
similar to the well-known viscous accretion discs, except 
for the changed sign of the rate of the mass flow. Angu- 
lar momentum is transferred from the optical compan- 
ion to the inner edge of the disc, increasing its angular 
velocity to Keplerian (rotationally supported). The ra- 
dial velocity component is subsonic even at the distance 
where the neutron star lies. Truncation occurs by the 
tidal interaction when the resonant torque exerted by 
the neutron star exceeds the viscous torque. This oc- 
curs only at certain radii — where the ratio between 
the angular frequency of disc rotation and the angu- 
lar frequency of the mean binary motion is a rational 
number. The efficiency of truncation depends strongly 
on the gap, Sr, between the truncation radius and the 
inner Lagrangian point and the viscosity parameter a. 
If Tdrift is the time scale for a particle in the disc to 
cross this gap then truncation will be efficient when 
Tdrift > ^orb- Typical values of the viscosity parameter 
are a < 1. Truncation is expected to be more efficient 
in low and moderate eccentricity systems (e<0.3) with 
narrow orbits (Porb ^ 40 d) than in high eccentricity 
sy stems because t he for mer have wider gaps. 



Okazaki et al.l (2002) performed numerical simula- 



tions of the disc formation around isolated Be stars and 
of the interaction between the Be-star disc and the neu- 
tron star in BeXB and showed that resonant truncation 
is effective if the viscosity parameter is a < < 1 . There 



is a radius outside of which the azimuthally-averaged 
surface density decreases steeply. When a ~ 0.1 the 
slope outside this radius is steeper, giving a stronger 
truncation effect on the disc. In the truncated disc 
model, the material lost from the Be star accumu- 
lates and the disc becomes denser more rapidly than 
if around an isolated Be star. 

5.2 Correlations with orbital period 

If the circumstellar disc is responsible for the variability 
observed at all wavelengths in BeXB, then correlations 
of an observational parameter with the orbital period 
can be attributed to the interaction of the neutron star 
with the disc. 

5.2.1 The Porh-EW(Ha) diagram 

Figure [15] shows an updated yersioii of the Porb~ 
EW(Ha) diagram by Reig et al.l ( 1997 . see also iReij 
(|2007t ): lAntoniou et al.l |2009bl )~. The original Porb" 
EW(Ha) diagram contained 11 BeXB, of which only 9 
had well established orbital periods. Since then sev- 
eral new BeXB have been discovered and their or- 
bital periods determined. Especially remarkable has 
been the discovery o f a large numbe r of BeXB in the 



Magellanic Clouds JCoe et al. 2005 : Liu et al 



Schmidtke fc Cowlvl l200a lAntoniou et al.ll2009bD . In 



2005 



this figure, only sources with wcU-detcrmincd orbital 
periods have been included, thus other SMC X-ray pul- 
sars (such as SXP756, SXP8.80, SXP46.6, and SXP304) 
with published Ha equivalent widths arc not present in 
the figure. 

The Porb~EW(Ha) diagram is based on the fact that 
the Ha line is the prime indicator of the circumstel- 
lar disc state. Although an instantaneous measure- 
ment of the EW(Ha) may not be an effective mea- 
surement of the size of the disc, the maximum equiv- 
alent width, when monitored during a long length of 
time (longer than the typical time scales for changes 
in the circumstellar disc, namely, from few months 
to a few years) becomes a significant indicator of the 
size of the disc. Based on int crferomctric observa tions, 
Quirrenbach et al.l (Il997l ) and lTvcner et al.l (|2005l) have 
shown that there is a clear linear correlation between 
the net Ha emission and the physical extent of the emit- 
ting region. Also, the observed correlations between the 
spectral parameters of the Ha line (FWHM, EW, peak 
separation in double-peak profiles) and the rotational 
velocity that have been observed in many Be stars are 
interpreted as e vidence for rotation ally dominated cir- 
cumstellar disc ( Dachs et al.lll986() . In particular, in- 
terpreting the peak separation (Apeak) of the Ha split 
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profiles as the outer radius ( Rput) of tlie emission line 
forming region (lHuang|ll97j) 



Rn 



2v sin i 



A. 



peak 



(1) 



the radius of the ernitting region can be e stimated 
(jHummel fc Vranckenlll995t IJaschek fc Jaschek.2004) - 
V sin i is the projected rotational velocity of the B star 
{v is the equatorial rotational velocity and i the in- 
clination toward the observer). As the EW(Ha) in- 
creases, the peak separation decreases, h ence iiicreas- 
ing the radius o f the Ha emitting region ( Dachsl[l992 : 
HanuschikI 119881) . Single peak profiles typically corre- 
spond to larger discs. 

Note that the value of the Ha equivalent width of 
Fig. [15] is the maximum ever reported. Some systems, 
especially in the SMC systems, have been discovered 
recently, hence the monitoring of the Ha is necessarily 
short. For these systems, the value of the EW(Ha) 
may not represent a high state of the circumstellar disc. 
The EW(Ha) of such systems should be considered as 
lower limits. This is the case of the galactic BeXB GRO 
J2058-h42 and most likely of the SMC source SXP455 
(RX JOlOl.3-7211). 

Assuming then, that EW(Ha) provides a good mea- 
sure of the size of the circumstellar disc. Fig. [15] in- 
dicates that systems with long orbital periods have 
larger discs, while narrow orbit systems contain smaller 
discs. The existence of the above relationship provides 
a strong argument that the discs are truncated, due to 
the action of the compact object, to a radius that is 
a fraction of the scmimajor axis. Given the similar- 
ity in masses of the optical companions in BeXB, we 
would expect that the truncation radius rtrunc would 
scale with the semimajor axis according to Kepler's 
Third Law as rtrunc oc P^/^. Since EW(Ha) scales 
as the square of the disc radius all other things being 
equal (jCrundstrom fc Gieslliooei) . the simple prediction 
would be EW{Ha) cx P'^^^. The solid line in Fig. [H] 
corresponds to a simple linear regression fit, while the 
dashed line represents a power-law fit with exponent 
equal to 4/3. 



5.2.2 Be 



B 



Be phase 



Column 4 of Table [5] gives the duration of the disc 
formation/dissipation phase, Tdisc, of those BeXB for 
which reliable data exist. The orbital period is given in 
column 6. As can be seen there is a good correlation be- 
tween Tdisc and the orbital period. Systems with short 
orbital periods tend to show faster disc growth and dis- 
sipation cycles, while longer time scales are associated 



with longer orbital periods. This result agrees with 
the scenario that the Be disc is truncated by the tidal 
torques from the neutron star, which would be more 
significant for narrower orbits. 

5.2.3 Long-term X-ray variability 



Reigj (|2007r ) investigated the long-term X-ray variabil- 
ity, expressed as the root-mean-square (rms) of the 
kSM/RXTE hght curves, of a set of BeXB and found 
that the systems with the larger rms are those har- 
bouring faster rotating neutron stars and those in low 
eccentric and relatively narrow orbits. The rms am- 
plitude was computed as rms = cr^/a;^, where x is the 
mean count rate and a"^ = a, 



obs 



j2 

cxp 



between the observed variance, 0-^^^^ = '^^{xi — x)^/N, 
and the expected variance, cr^^p = '^i^ll^ {^i are 
the experimental errors, and TV is the total number of 
points). 

Figure [16] shows the rms as a function of the spin 
period, eccentricity and orbital period of the systems. 
Note that by using the rms the emphasis is put on the 
amplitude of the variations rather than on how often 
the variations take place. For the sake of comparison, 
persistent BeXB and supergiant X-ray binaries, that is, 
systems whose optical companion is an evolved (lumi- 
nosity class TII) star, have also been included. Different 
type of systems have been represented by different sym- 
bols as follows: transient BeXB by circles, persistent 
BeXB by stars and supergiants by squares. Systems 
that have shown Type II outbursts are represented by 
black filled circles. Triangles denote the class of low-e 
BeXB (Sect. \Un\) . It is very likely that low-e BcXB 
and persistent BeXB form only one class of systems. 
In fact, X-Per, which is considered the prototype of 
the persistent BeXB, also belongs to the class of low-e 
BeXB. 

In all cases a clear anticorrelation is apparent: the 
systems with the faster rotating neutron stars and with 
low eccentricity and narrow orbits are more variable, 
i.e., present higher rms. This plot then suggests two in- 
teresting results. First, since the BeXB and the SGXB 
occupy clearly distinct regions in Fig. 1161 the X-ray 
variability in high-mass X-ray binaries does not only 
depend on the physical conditions in the vicinity of the 
compact object but also on the mass transfer mech- 
anism, i.e., whether a circumstellar disc is present or 
not. Second, the systems containing fast spinning neu- 
tron stars are more likely to exhibit Type II outbursts. 

In general, i) the more eccentric the orbit the lower 
the variability and ii) type II activity mainly occurs 
in low-eccentric systems. The low eccentricity class of 
BeXB follow the general trend in the rms — Pspin and 
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Fig. 16 X-ray rms variability of tlie RXTE/ ASM light curves as a function of the spin period (left), eccentricity (middle) 
and orbital period (right) of galactic BeXB. Black filled circles denote systems that have shown type II outbursts. Squares 
represent SGXBs and the star symbols corr espond to pe rsistent BeXB. The cyan filled circle on the right panel correspond 
to SAX J2103. 5+4545 (see Sect.[ri|). From lReid (|2007l ). 



rms— Porb plots but clearly distinguish themselves from 
the classical BeXB, in the rms — e diagram. The an- 
ticorrelation of the rms with the orbital period pro- 
vides further evidence for the truncation of the disc in 
systems with short orbital periods, i.e., those having 
narrow orbits. 

A prediction of the viscous dccretion disc model is 
that type I outbursts are expected to occur more of- 
ten in high-eccentric systems, whereas type II outburst 
would dominate in low-eccentric systems, in agreement 
wi th the resu l ts sho wn in Fig. [T51 



Coe et all (120051 ) carried out a major study of the 



optical and infrared characteristics of BeXB in the 
SMC. They found that the optical photometric vari- 
ability is greatest when the circumstellar disc size is 
least and also in syste ms containing fast rotating neu- 



tron stars. iReia (|2007[ ) extende d these results to the X- 
ray do r nain. While the work s bv lOkazaki fc Negueruela 
I 200l[ ). IZhang et all (|2004[ ) and the EW{Ha) - P^b 
correlation clearly identify truncation with low eccen- 
tric and narrow orbit systems, the plots of the X-ray 
rms as a function of the system parameters (Fig. 1161) 
and those of the EW(IIq;) and Pspin as a f unction of 
the o ptical photometric rms (Figs 6 and 7 in ICoe et aF 



20051) identify truncation with variability. That is, sys- 
tems in which truncation is favoured are the most vari- 
able ones in both, the X-ray and optical bands. 

5.3 Be in BeXB versus isolated Be 

The other line of argument to demonstrate that the 
neutron star exerts a measurable effect on the optical 



companion is based on the comparison of the variability 
patterns of the optical/IR emission in isolated Be stars 
and BeXB. 

5.3.1 V/R variability 

As explained in Sect. 12.11 V/R variability is defined as 
the intensity variations of the two peaks (known as vi- 
olet and red peak) in the split profile of a spectral line. 
In many Be stars, if monitored over a long en ough pe- 
riod o f time, these variations are quasiperiodic (jOkazaki 
19971 ). V/R variability is currently explained by the 
global one-armed oscillation model. This model sug- 
gests that the long-term V/R variations are caused by 
global m = 1 oscillations in the cool equatorial disc of 
the Be star. In other words, an enhanced density per- 
turbation develops on one side of the disc, which slowly 
precesses. The precession time is that associated with 
the V/R quasi-period. The periods of V/R variations 
in isolated Be stars range from years to decades, with a 
statistical mean of 7 years, which is m uch longer tha n 
the rotational period of the central star ( Okazakilll997 ). 
In contrast, in BeXB (Table [5]), the V/R quasiperiods, 
Ty/Rj are significantly shorter (<5 yr ), which can be 
interpreted as the result of smaller discs in BcXB due 
to truncation caused by the ne utron star. 

Oktariani fc Okazaki (I2OO9I) studied the tidal effect 



of the companion on the global oscillation modes in 
equatorial discs around binary Be stars. They found 
that the oscillation modes are well confined when the 
disc is larger than a few tens of stellar radii. In smaller 
discs, however, the mode confinement is incomplete and 
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Fig. 17 Peak separation as a function of EW(Ha). Tlie 
daslied line represents tiie average beiiaviour of tiie iso- 
lated Be stars. Be/X-ray binaries show, on average denser 
disks. Different sysmbols correspond to different BeXB. 
From lZamanov et al.l (|200fl '). 



the oscillation period depends on the binary parameters 
significantly. The oscillation period becomes longer in 
the case of a wider binary separation and/or a lower 
binary mass ratio. Such dependence is observed only 
for short period binaries. In wide binaries where the 
disc size is large enough to confine the one-armed modes 
in the inner part of the disc, the mode characteristics 
depend little on the binary parameters. 

5.3.2 Disc density 



Zamanov et al.l (|2001l ) carried out a comparative study 
of the circumstellar discs in BeXB and isolated Be stars 
based on the Ha line emission properties. They found 
that the discs of BcXB are twice as dense as and have, 
on average, smaller size than the discs of isolated Be 
stars. 



HanuschikI ( 19881 Il989f ) showed that there is a corre- 



lation between the Ha peak separation and the equiv- 
alent width in Be stars, according to the law 
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where usini is the projected rotational velocity and 
EW(Ha) is given in Angstroms, a and b are related to 
the rotational law index a = — j/2 (j=0.5 for Keplerian 
rotation and j=l for conservation of angular momen- 
tum) and with the disc electron density, respectively. 
Figure [JT] shows the plot of log ^^Yini '"^ EW(Ha) 
for a number of BeXB. The straight line (a = —0.4, 
b — —0.1) represents the aver age behavio u r of t he iso- 



lated Be stars investigated by iHanuschikl (|l989l ). The 
data points correspond to BcXB. As can be seen, the 



vast majority of the data points lie above the average 
line, i.e. they are shifted towards denser discs. 

Another result that supports denser disc in Be bina- 
ries is the correlation observed between the Ha equiv- 
alent width and the (B — V) colour excess in a sam- 
ple of Be/X-ray bin aries in the Magellanic Clouds 
( Torreion et al. 2007t ). The total reddening in the di- 
rection of a Be star or BeXB binary is composed of 
two components: one (constant) produced mainly by 
dust in the interstellar space through the line of sight 
£'^"{3 — V) and another one (variable) produced by 
the circumstellar gas around the Be star E^^{B — V). 
The nature and wavelength dependence of these two 
reddenings is completely different but their final effect 
upon the photometric ind ices and spectrum is, in prin- 



ciple, indistinguishable. iTorreion et al.l (120071 ) found 



that the circumstellar excess is ~ 10 times larger for 
Be stars in binaries than for isolated Be stars with the 
same amount of emission, and also that the circumstel- 
lar colour excess, E^^{B — V) saturates at~ 0.2 magni- 
tudes beyond EW(Ha)«20 A. 

Again the denser discs seen in BeXB can be ex - 



plained by disc truncation. lOkazaki et al.l (|2002l ) 
showed that during the initial build-up phase, the disc 
evolution is similar to that for isolated Be stars. But, 
later on, the effect of the resonant torque becomes ap- 
parent, preventing the disc gas from drifting outwards 
at several resonance radii. As a result, the disc density 
increases more rapidly than that for isolated Be stars, 
in agreement with the results presented in this section. 

5.3.3 Disc size 

As mentioned above, interferometric observations have 
shown that there is a clear linear correlation between 
the net Ha emi ssion and the physical extent of the 



emitt ing region ([Quirrenbach et al.lll997l : iTvcner et al 
20051) . Using nur rierical models of th e circu mstellar 



discs of Be stars, Grundstrom fc Gie^ ( 20061 ) showed 
that there are monotonic relationships between the 
emission-line equivalent width and the ratio of the an- 
gular half-width at half maximum of the projected disc 
major axis to the radius of the star, providing a method 
to estimate disc radii based only on the Ha equivalent 
width with an accuracy of about 30%. Since, on av- 
erage, isolated Be stars show larger EW(Ha) than Be 
stars in BeXB |Reig et al.ll997tlAntoniou et al.i2009bl) , 
we conclude that, on average, the envelopes of Be stars 
in BeXB are smaller than those present in isolated Be 
stars. 

5.4 Quantized infrared excess flux states 

Further evidence in favour of the decretion disc model 
and disc truncation comes from the infrared band. 
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Table 6 Characteristic timescales for disc variability in Be/X-ray binaries. 



X-ray 


Optical 


Spectral 


^disc 




-Porb 


Reference 


name 


name 


type 


(year) 


(year) 


(day) 




4U 0115+63 


V635 Gas 


B0.2Ve 


3-5 


0.5-1.5 


24.3 


1 


V 0332+53 


BQ Cam 


09Ve 


4-5 


1 


34.2 


2,3 


4U 0352+309 


X Per 


BOVe 


7 


0.6-2 


250 


4 


lA 0535+262 


V725 Tau 


09.7IIIe 


4-5 


1-1.5 


111 


5,6 


RX J0812.4-3114 


LS 992 


B0.2IVe 


4 




80 


7 


SAX J2103.5+4545 




B0.5V 


1.5-2 




12.7 


8 


RX J 0146.9+6121 


LS I +61 235 


BlVe 


>10 


3.4 


>200* 


9 


4U 1145-619 


V801 Ccn 


BlVe 


>10 


3 


186 


10 


RX J0440.9+4431 


LS V +44 17 


B0.2Ve 


>10 




>150* 


11 


IGR J06074+2205 




B0.5V 


>5 


-5 




12 



* Obtained from the Pg, 



correlation. 



Note: fl] [Negueruela et al.l 1I2OO1I'). [2] iNegueruela et al l \l99^ ). f3 |lGoranskii[ l [200l|) f4] |Clark et all 1I2OOII). [S l lHaigh et al.l \2004). [6 ] 
ICIark et al.l lll998l') [7] Reig et alj l200lh . fSl lReig et al.l l l2010af ). [Ql lReig et ahl i200ol ) flO] IStevens et al.l lll997l '). [11] iReig et al.l ||2005| '). 
[12] iReig et al.l ll2010bli 



Haigh et al. ( 20041 ) showed that the Be/X-ray binary 
A 0535+26 displays a clear bimodality in the infrared 
magnitudes distribution, suggesting that the system al- 
ternates between a faint and a bright flux state. 

These quantized infrared excess flux states are a nat- 
ural consequ ence of the resonant tru r rcatiq n hypothesis 
proposed bv lOkazaki fc Negueruelal (|2001l ) and consti- 
tute a firm observational verification of the proposed 
resonant truncation scheme for Be X-ray binaries. Ac- 
cording to this model, tidal torques, which operate most 
strongly in orbits resonant with that of the neutron 
star, truncate the disc at radii meeting t he condition 
Pns ■ Ptrunc ~ u : 1 whcrc n is an integer. iHaigh et al. 
( 20041 ) identified the two states with resonances 5:1 and 
6:1. A possible brighter infrared state corresponding to 
resonance 4:1 may have also been seen during the 70's. 

The simultaneity of X-ray activity with transitions 
between these states strongly suggests a broad mecha- 
nism for outbursts, in which material lost from the disc 
during the reduction of the truncation radius is accreted 
by the neutroir star. 

The idea is simple and neat. When the disc outer ra- 
dius makes a quantum leap inwards the material present 
between those radii must be relocated elsewhere in the 
system and becomes available for accretion. In 1993/94 
a change of infrared state as accompanied by a giant 
X-ray outburst. This simultaneity of the type II out- 
burst with the change of infrared state strongly suggests 
that the giant outburst was caused by the accretion of 
material previously resident between the 5:1 and 6:1 
resonances. Importantly, truncation would appear to 
be the cause of the X-radiation and not vice versa, as 
truncation commenced before X-ray emission. 



6 Be/X in other galsixies 

The study of X-ray source populatioirs in other galaxies 
is hampered by the large distance from us. Even other 
galaxies of the Local Group are too far to reach the 
quiescence population of HMXBs (Lx ^ 10^^ — 10^^^ 
erg s ~^). Good enough signal-to-noise optical spec- 
tra, needed to derive the spectral type of the optical 
counterparts, are very hard to achieve. Because of the 
relatively low intergalactic extinction, well-established 
distance and their proximity, the Magellanic Clouds 
(]VIC) represent the best targets to study a complete 
X-ray binary population. In fact, they are the only 
galaxies for which meaningful statistical studies can be 
performed (current space detectors can detect MC X- 
ray sources down to a few 10'^'^ erg s ^^). Moreover, the 
MC show diffcrcirt chemical compositions and are heav- 
ily interacting with the Milky Way, which affects their 
star formation history. Therefore, the study of stellar 
populations in the MC is particularly rewarding. 

Although past X-ray i nissions with imag ing capa- 
bilities, such as Ei nstein dWang fc Wulll992l). ROSAT 



( Haberl et aflboOOl) . and ASCA (|Yokogawa et allliooi) 



have surveyed the MC, the deep surveys of the current 



20081). CHANDRA (I 


( Galache et al. 


20081) 


20061: Coe et al. 


2010l) 



(|McGowan et al.l 120081). RXTE 



and INTEGRAL ( Gotz et al 



complemented with observa- 



band (Schmidtke & Cowlv 2006 


; Schurch et al. 


2007; 


McBride et al. 20081 Antoniou et al. 2009aibl) are the 



ones giving new insights into population synthesis mod- 
els. 
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The Small Magellanic Cloud (SMC) hosts an unex- 
pected large population of HMXBs. Initial estimates 
of the number of HMXBs in the SMC were based on 
the mass of the SMC relative to that of the Milky Way. 
The SMC is a few percent of the mass of the Galaxy, 
in which about 70 X-ray pulsars are known. Therefore 
one would expect to find only two or three systems. 
However, there are about 50 HMXBs in the SMC. Cu- 
riously, only one (SMC X-1) of these systems is not a 
Bc/X, that is, it contains a supergiant companion. 

The reduced metallicity of the SMC one fifth of 
solar) cannot explain, by itself, t he large number of 
HMXBs in this galaxy (jPravlbood) . In addition to re- 
duced metallicity, a large star forr nation rate is neede d 
to explain this excess of HMXBs ([Grimm et al.l 120031 ) . 
The star formation rate/mass of the SMC is 150 times 
that of the Milky Way. The star formation rate may 
have been enhanced as a result of tidal interactions be- 
tween the MC and the Galaxy. However, the most re- 
cent close approach of the SMC and LMC was 200 Myr 
ago - much longer than the evolutionary timescale of 
Be/X-ray binaries. This implies that either there has 
been a significant delay between the encounter of the 
SMC and LMC and the onset of star formation, or that 
subsequent waves of star formation have given r ise to 
these Be/X-ray binaries ( Harris &: Zaritskvl 20041 ) . 

The main question that we wish to address in this 
section is whether the population of BeXB in the MC 
differs from that of the Milky Way. We focus on the 
statistical distributions of spin periods, spectral type 
and Ha equivalent width of the optical counterparts. 

6.1 Spectral-type 

The spectral distribution of the BeXB population in 
the Milky Way is characterised by a sharp peak at 
BO and a relatively narrow width, namely, no BeXB 
are found beyond type B2 and earlier than type 
08. iNegueruela fc Coel (|2002[) : lAntoniou et al.l |2009bl) 
showed that BeXB in the LMC follow the same dis- 



tribution as those in our Ga l axy, w hile iMcBride et al 



20081) and lAntoniou et al.l (|2009bl ) showed that the 



spectral distribution of BeXB in the SMC is also con- 
sistent with that of the Milky Way. The only difference 
might be the presence of a very small number of BeXB 
with spectral type later than B2 in the SMC. However, 
given the very low number of these systems and the 
uncertainty in the spectral classification, this difference 
is not statistically significant. 

6.2 Ha equivalent width 

On average, the Ha equivalent width (EW(Ha)) of the 
SMC population of BeXB is larger than that of the 



-- SMC 

— Milky Way 



10 - 



I 5- 



-1 1 2 3 

l°g(Ppi„) 

Fig. 18 Distribution of pulse periods in the SMC and the 
Galaxy. 



Milky Way. lAntoniou et al.l (|2009bl ) found that the 
EW(Ha) distribution of SMC BeXB peaks in the in- 
terval from -15 to -25 A, while the EW(Ha) of Galac- 
tic BeXB distribute more uniformly in the interval 
from -30 to A. Likewise, the number of sources with 
EW(Ha) > 50 A is larger in the SMC. For the LMC sys- 
tems the sample is too small to draw any conclusions. 
This trend may just be a selection effect. The discovery 
of BeXB in the SMC results mainly from X-ray obser- 
vations. Due to their larger distance we tend to detect 
the brightest systems, hence those with a large reservoir 
of accreting material, i.e., those with larger discs. 

6.3 Spin periods 

As indicated in Sect. [21 virtually all HMXBs are X-ray 
pulsars. Histograms of the spin period of HMXBs in 
the MC and the ir comparison with G alact ic sources 
can b e found in Lavcock et al. ( 20051 ) and Liu et al 
(|2005l ). They both find that the SMC population ap- 



2005) 



20051 ) 



pears to be shifted to shorter pe riods. However , while 
this shift is very significant in the iLavcock et al 
sample, the difference is minor in the lLiu et al. 
sample. One contribution to the difference between 
the Galactic and SMC pulse period distributions is the 
apparent lack of supergiant wind accretion pulsars in 
the SMC. HMXBs with supergiant companions tend 
to have longer (> 100 s) spin periods than BeXB. As 
stated before only one such system is known in the 
SMC. Figure [T8| shows the spin period distributions 
of well estabhshed BeXB (i.e. excluding supergiant 
wind-powered systems) in the SMC and Milky Way 



from the lists comp i led b y iLavcock et al.l ([20051 ) and 
Raguzova fc Popovl (|2005l ). A Kolmogorov-Smirnov 



test gives no real evidence against the null hypothesis 
that the two samples are drawn from the same distri- 
bution. 
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7 Summary 

The traditional picture of two classes of high-mass 
X-ray binaries, namely supergiant X-ray binaries and 
Be/X-ray binaries is giving way to a more complex 
situation where newly discovered systems may not fit 
in these categories. Superfast X-ray transients, low 
eccentricity BeXB, 7-ray binaries, and 7-Cas like ob- 
jects are among the new type of systems that only re- 
cently have begun to emerge. BeXB do not form a ho- 
mogeneous group either with transient and persistent 
sources, highly eccentric and nearly circular orbits, and 
fast (Pspin ^ few seconds) and slow {Pspin ^ few hun- 
dred of seconds) rotating neutron stars. 

Virtually, all BeXB are X-ray pulsars. The study 
of the variability of the X-ray spectral and pulse tim- 
ing parameters across giant outbursts is crucial to un- 
derstand the physics of accretion in strongly magne- 
tised neutron stars. Only recently, attention is being 
paid to the aperiodic variability and evolution of broad- 
band noise components across the outburst. The ap- 
plication of diagrams such as hardness-intensity and 
colour-colour diagrams, so widely used in black-hole 
and low-mass binaries, to BeXB and the definition of 
spectral states is a promising way to understand the 
phenomenology of the evolution of outbursts in these 
systems. 

Progress in understanding BeXB is also being made 
in the optical and infrared bands where the massive 
companion shines bright. The data gathered over the 
past ten years indicate that the interaction between the 
compact object and the Be type star works in two di- 
rections: the massive companion provides the source 
of matter for accretion; the amount of matter captured 
and the way it is captured (transfer mechanism) change 
the physical conditions of the neutron star (e.g. by 
spinning it up or down and shrinking or expanding the 
magnetosphere). But also, the continuous revolution of 
the neutron star around the optical counterpart pro- 
duces observable effects, the most important of which 
is the truncation of the Be star's equatorial disc. 

Future missions with improved sensitivities in the X- 
ray band will help solve remaining open questions such 
as the formation of accretion discs during X-ray out- 
bursts, the origin of quasi-periodic oscillations, charac- 
terisation of the population of persistent BeXB and the 
triggering mechanism of the giant outbursts. 
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